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PERMITTIVITY  MEASUREMENTS  UTILIZING  THE 
DIELECTRIC  POST  IN  WAVEGUIDE  TECHNIQUE 

SUMMARY 

This  report  describes  experimental  procedures  and  computer  codes 
developed  at  NRL  for  the  measurement  of  the  permittivity  of 
filamentary  high  loss  materials  by  the  dielectric  post-in-waveguide 
(DPW)  technique.  This  method  is  especially  useful  in  characterizing 
the  microwave/millimeter  wave  dielectric  properties  of  materials  such 
as  partially  graphitized  carbon  fibers,  linear  chain  crystals, 
superionic  conductors,  ferroelectrics,  etc.  This  method  allows  the 
complex  permittivity,  to  be  measured  on  a  pointwise  continuous 
frequency  basis  across  a  rectangular  waveguide  band.  We  have 
employed  and  will  describe  two  measurement  setups;  for  X  and  Ku  band 
studies  we  have  collected  data  with  commercial  automated  vector 
network  analyzers,  for  K  and  Ka  bands  we  have  assembled  transmission 
phase  bridges  using  standard  waveguide  components  and  a  scalar 
network  analyzer.  The  data  from  both  measurement  configurations  are 
collected  on  floppy  disks  for  later  analysis  on  a  VAX  11/780 
computer.  The  computer  codes  developed  for  this  work  will  be 
described  in  detail  and  selected  programs  and  subroutines  will  be 
listed  in  the  appendix. 


INTRODUCTION 

There  are  many  methods  of  measuring  the  dielectric  constant  of 
materials  at  microwave  frequencies.  The  simpliest  and  the  most 
broadband  technique  employs  an  automated  vector  network  analyser 
( AVNA )  and  a  coaxial  or  waveguide  transmission  line.  The  sample  is 
chosen  to  be  somewhat  less  than  a  quarter  wavelength  thick  at  the 
lowest  measurement  frequency  and  to  fill  the  cross  section  of  the 
transmission  line.  The  AVNA  measures  the  S-parameters  of  the  loaded 
transmission  line  and  the  frequency  dependent  permittivity,  e(f),  and 
permeability,  u(f),  of  the  sample  can  be  directly  calculated  from  the 
S-parameters  by  the  AVNA  computer/controller.  This  technique, 
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however,  assumes  that  the  sample  is  homogeneous  in  composition  and 
isotropic  in  microwave  dielectric  properties.  These  assumptions  fail 
for  highly  anisotropic  samples.  At  best,  Mie  filled  transmission 
line  technique  can  be  employed  for  measuring  anisotropic  materials  in 
the  special  case  where  the  anisotropy  axes  of  the  pure  material  can 
be  aligned  with  the  internal  fields  of  the  transmission  line  or  to 
measure  the  composite  properties  of  anisotropic  materials 
mechanically  dispersed  in  a  binder.  In  the  case  of  the  composite 
system  it  is  very  difficult  to  infer  the  intrinsic  properties  of  the 
anisotropic  material  from  the  composite  dielectric  response. 

A  second  technique,  cavity  perturbation,  is  a  classic  method  for 
measuring  the  permittivity  and  permeability  of  materials  and  is 
suited  for  the  study  of  filamentary  samples.  The  perturbation 
technique  ’fails  for  high  loss  samples  which  destroy  the  cavity 
resonance.  Reduction  in  sample  size  can  alleviate  the  problem  in 
some  cases,  however,  impossibly  small  and  very  fragile  samples  with 
large  geometric  aspect  ratios  may  be  required  to  satisfy  the 
perturbation  criterion  while  minimizing  depolarization  effects. 
Perturbation  cavities  also  have  the  disadvantage  of  operating  at  a 
single  frequency. 

The  dielectric  post-in-waveguide  method  for  measuring  the 
permittivity  of  filamentary  or  cylindrical  post  samples  shares  some 
of  the  advantages  of  the  two  conventional  methods  previously 
mentioned.  Like  the  filled  transmission  line  method  it  allows 
relatively  broadband  automated  measurments  and  like  the  perturbation 
cavity  is  it  is  adapted  for  the  determination  of  permittivities  of 
filamentary  samples.  The  technique  works  best  with  samples  with 
large  dielectric  constants  and/or  high  loss  -  samples  which  fail  in 
the  resona't  cavity  perturbation  technique.  Samples  chosen  for  study 
by  this  technique  must  present  a  sufficient  impedance  mismatch  to  the 
waveguide  to  measurably  affect  the  reflection  and  transmission 
coefficients.  Since  in  this  method  the  cylindrical  filament  axis  of 
the  sample  is  alligned  with  the  microwave  E-field,  the  permittivity 
measured  is  that  of  the  long  axis  of  the  sample. 

The  method  waa  first  described  by  J.  Le  Bot  who  applied  it  to  the 
measurement  of  e  as  a  function  of  water  absorption  in  silica  and 
alumina  gels  (1).  Since  that  time  the  method  has  been  used 
relatively  infrequently  (2,3).  In  this  report  we  describe  how  this 
technique  can  ~acome  a  useful  tool  for  research  using  modern 
instrumentation  and  computational  capabilities. 

The  QPW  experimental  procedure  is  greatly  simiplified  by  the  use 
of  the  automated  vector  network  analyser.  AVNA's  are  now 
commercially  available  for  pointwise  continuous  swept  frequency 
measurements  up  to  frequencies  of  100  GHz.  The  AVNA  available  to  our 
group,  however,  is  limited  to  frequencies  less  than  or  equal  to  18 
GHz.  He  have  thus  employed  the  AVNA  for  post  measurements  at  X  and 


Ku  bands.  Measurements  were  extended  to  K  band  (18  0  26.5  GHz)  Ka 
band  (26.5  to  40  GHz)  by  the  construction  of  simple  transmission 
bridges.  In  either  case  the  measurement  fixture  consist  of  a  simple 
short  length  of  waveguide  with  a  hole  drilled  thru  the  center  of  the 
broad  face.  (See  Fig.  1).  The  waveguide  sample  holder  can  be  heated 
in  a  tube  furnace  or  cooled  to  about  100K  by  standard  gas  flow 
techniques.  The  swept  frequency  data  is  collected  on  floppy  disks  on 
both  the  AVNA  and  transmission  bridge  set-ups  and  analyzed  on  a  VAX 
11/780  computer. 

The  computer  codes  created  for  this  analysis  employ  subroutines 
which  calculate  the  complex  transmission  and  reflection  coefficients 
for  the  post  geometry.  One  code  calculates  the  return/insertion  loss 
and  phase  expected  for  a  given  frequency  dependent  permittivity, 
samde  diameter,  and  waveguide  width.  Another  program  uses  the 
measured  return  loss  and  phase,  or  insertion  loss  and  phase,  or 
scalar  return  loss  and  insertion  loss  to  calculate  permittivity  of  a 
given  sample  at  each  measurement  frequency  across  the  waveguide  band 
using  a  Newton  inversion  technique.  A  third  code  performs  a  least 
square  fit  to  S-parameter  data  at  each  frequency.  A  fourth  code 
performs  at  least  squares  fit  to  multiple  sets  of  S-parameter  data 
assuming  a  menu  determined  frequency  dependent  functional 
representation  of  the  complex  permittivity  of  the  sample  over  one  or 
more  waveguide  bands. 

In  the  following  sections  we  discuss  the  theory,  the  experimental 
set-ups,  the  computer  codes,  and  provide  some  examples  of  the 
measurement  technique.  In  the  appendices  we  provide  a  detailed  list 
of  the  measurement  and  analysis  procedures  and  list  some  of  the 
computer  codes  developed  for  DPW  data  reduction. 


THEORY 

The  problem  of  a  dielectric  post  in  a  rectangular  waveguide  was 
solved  by  workers  of  the  MIT  Radiation  Lab  during  World  War  II.  The 
method  of  solution  is  discussed  in  detail  in  a  volume  of  lecture 

notes  from  that  era  by  Prof.  Julian  Schwinger  (4)  and  the  results  of 
the  calculation  to  high  order  were  presented  in  the  MIT  Waveguide 
Handbook  (5).'  In  most  experimental  studies  employing  the  post  method 
it  is  assumed  that  the  sample  diameter  is  very  small  with  respect  to 
the  microwave  wavelength  and  waveguide  dimensions.  In  this  limit  the 
equations  can  be  simplified  somewhat  and  permittivity  may  be  found 
using  graphical  methods  (1).  Here  we  assume  only  that  the  sample  is 
cylindrical  as  required  by  the  theory  and  that  it  is  positioned  in 
the  center  of  the  broad  face  of  the  waveguide  with  its  cylindrical 

axis  parallel  to  the  microwave  E-field. 

The  dielectric  post  problem  is  analyzed  in  terms  of  a  T  shaped 

equivalent  circuit  where  the  two  in  line  impedance  elements,  Zb,  on 
top  of  the  T  are  identical  and  differ  from  the  shunting  impedance, 
Za.  The  circuit  is  assumed  to  be  terminated  on  each  end  by  the 
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characteristic  impedance  of  the  waveguide  Zo*  (See  Fig.  2). 
this  case  the  equations  become: 
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where  a  is  the  waveguide  width  and  d  is  the  sample  diameter  A  is  the 
free  space  wavelength  and  is  the  wavelength  in  waveguide.  The 
computer  calculation  uses  the  full  expression  since  some  samples  with 
large  dielectric  constants  are  not  electrically  thin.  The  sum  is 
taken  to  100  terms  for  odd  n.  The  complex  J0  amd  Ji  Bessel 
functions  are  evaluated  by  series  approximations  (6).  The 
S-parameters  are  given  by: 

sn  -  s22  =  [(z1  +  z2)2  -Z2  -i]/[(z2  +  zx  +  l)2  -  z|]. 


and 


Si 2  -  S21  =  2Z2/[Z?+2Z2Zi+2(Z2+Zi)+1]., 

where  standard  waveguide  theory  was  used  to  derive  the  measured 
circuit  parameters  from  the  equivalent  circuit  representation  (7). 


EXPERIMENTAL  SET-UPS 

DPW  measurements  have  been  greatly  simplified  with  the  advent  of 
accurate  automated  vector  network  analyzers.  These  devices  can 
rapidly  scan  the  waveguide  band  measuring  the  S-parameters  of  the  DPW 
waveguide  fixture  at  hundreds  of  frequencies.  The  X ( 7 . 5  to  12.4  GHz) 
and  Ku  (11  to  18  GHz)  data  were  taken  with  a  Hewlett-Packard  Model 
8409S  and  recently  with  a  Model  8510  automatic  vector  network 
analyzers.  For  measurements  with  the  8409S  AVNA  we  find  it 

convenient  to  terminate  the  waveguide  test  fixture  with  a  good 
matched  load  and  employ  the  Sll  parameter  for  the  computer 
inversion.  Using  just  the  Sll  values  speeds  up  the  measurement  time, 
especially  saving  time  in  the  initial  calibration  of  the  AVNA.  We 

find  it  necessary  to  place  a  pad  (about  5  dB)  before  the  waveguide 

fixture  to  damp  out  spurious  standing  waves.  Nevertheless,  the 

principal  source  of  error  arises  from  uncorrected  standing  waves  in 
the  measurement  arm.  With  the  acquisition  of  the  8510  AVNA  the  post 


Figure  2.  (a)  A  sketch  of  the  dielectric  post  positioned  1 
waveguide  where  a  Is  the  width  of  the  waveguide  and  d  Is  the  samp! 
diameter;  and,  (b)  a  sketch  of  the  equivalent  circuit. 
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attenuators  are  located  In  both  arras  of  the  bridge  with  a  calibrated 
phase  shifter  in  the  lower  (reference)  arm. 


measurements  are  much  improved  owing  to  an  order  of  magnitude 
improvement  in  the  accuracy  of  the  complex  S-parameter 

determinations.  With  the  8510  we  measure  both  Sn  and  S21  with 
small  standing  wave  distortions  employing  simple  unpadded  waveguide 
fixtures.  The  permittivity  can  then  be  infered  from  Sn,  S21  or 
Sn  cLn.d  S21  .  The  choice  of  data  type  is  determined  by  tht 
relative  reflection  and  transmission  amplitudes.  Least  squares 
fitting  of  the  data  employes  both  Sn  and  S21 . 

We  also  have  made  measurements  at  K ( 18  -  26.5  GHz)  and  Ka  band 
(26.5-40  GHz).  Our  AVNAs  do  not  cover  these  waveguide  bands  and  so 
we  have  assembled  simple  transmission  bridge  circuits  at  these 
frequencies.  (See  Fig.  3).  The  source  in  this  case  is  a  HP  Model 

8690B  Microwave  Sweeper  with  appropriate  BWO  plug-ins.  An  amplitude 
reference  signal  is  tapped  off  just  after  the  sweeper  and  then  the 
microwave  power  is  split  into  a  sample  arm  and  a  reference  arm.  Both 
sample  and  reference  arms  contain  precision  calibrated  attenuators 
and  the  sample  arm  holds  the  DPW  test  fixture  whereas  the  reference 
arm  contains  a  calibrated  phase  shifter.  Microwaves  from  the  two 
arms  are  brought  together  and  the  mixed  signal  is  detected  by  a 
standard  wave  guide  Shotkey  detector.  Measurements  are  acquired  by  a 
Model  1038-N10  Pacific  Measurements  Scalar  Network  Analyzer  and  fed 
to  a  LSI-11  computer  for  storage  on  floppy  disk.  Data  is  taken  as 
the  sweeper  is  swept  across  the  full  waveguide  band.  The  insertion 
loss  is  first  measured  with  the  reference  arm  sec  at  maximum 
attenuation.  Phase  data  is  then  acquired  with  the  sample  in  place 
and  with  an  empty  test  fixture.  For  the  phase  measurement  the  power 
in  the  two  arms  is  adjusted  so  that  the  power  through  the  reference 
arm  is  slightly  greater  than  that  through  the  sample  arm.  The  phase 
shifter  is  then  stepped  through  360  degrees  at  (usually)  20  degree 
intervals.  At  each  phase  setting  a  full  sweep  is  acquired  and  stored 
ori  the  LSI/11.  The  sample  arm  attenuator  is  usually  set  at  7  dB  to 
damp  out  spurious  standing  waves  and  an  isolator  is  used  before  the 
mixer-detector  for  the  same  purpose.  The  phase  data  is  then  analyzed 
with  a  computer  code,  PHASE,  whicn  performs  at  least  squares  fit  to 
data  at  all  the  phase  shifter  settings  at  each  acquisition  frequency 
(or  a  specified  subset  of  frequencies).  This  fitted  phase  shift  is 
normalized  to  the  fitted  phase  shift  for  data  taken  with  an  empty 
sample  holder.  In  some  cases  return  loss  is  also  measured  at  a 
second  detector  mounted  on  a  directional  coupler  in  the  sample  arm 
and  referenced  to  a  short  inserted  just  before  the  waveguide  test 
fixture. 

As  mentioned  earlier  the  waveguide  fixture  consists  of  a  simple 
short  length  of  waveguide  with  the  sample  inserted  through  a  hole 
centered  in  the  broad  face.  The  sample  is  usually  housed  in  a 
thin-walled  glass  capillary  tube  commonly  employed  in  Debye-Scherer 
X-ray  measurements.  Insertion  of  the  test  sample  by  means  of  a 
capillary  tube  through  the  hole  in  the  wave  guide  allows  rapid  sample 
insertion  and  for  a  minimum  disturbance  of  the  test  network.  This 
sample  mounting  also  allows  the  imposition  of  a  D.C.  voltage  across 
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the  sample  while  it  is  being  measured  at  microwave  frequencies.  This 
sample  insertion  method  may  cause  problems  in  the  case  of  high  or 
moderately  conducting  samples  or  for  high  dielectric  samples  due  to 
deoolan zation  effects  at  the  waveguide  walls.  For  such  samples  it 
is  probably  necessary  to  make  good  electrical  contact  to  the 
waveguide.  One  method  for  doing  this  is  to  use  a  waveguide  cover 
flange  with  a  small  grove  in  one  face  to  hold  the  sample.  The  sample 
is  cemented  in  place  in  the  grove  with  conducting  paste  (3).  On  a 
few  occasions  we  have  used  a  split  waveguide  fixture  to  hold  larger 
samples.  The  samples  are  cut  to  form  right  circular  cylinders  of 
very  nearly  the  waveguide  height  and  are  then  press-fitted  into  the 
split  section  of  waveguide  for  measurement. 


COMPUTER  CODES 

A  number  of  computer  codes  have  been  developed  to  simulate  and 
analyze  the  data  generated  by  the  dielectric  post  in  waveguide 
experiments.  The  major  codes  and  their  functions  are  described 
below.  Since  the  experimental  setups  are  simple,  most  of  the  effort 
in  creating  the  measurement  capability  was  devoted  to  developing  the 
computer  codes.  Likewise,  the  data  acquisition  is  simpler  than  the 
data  analysis,  especially  when  the  data  is  taken  with  the  automated 
vector  network  analyzer.  All  the  codes  described  in  this  report  are 
written  in  Fortran  IV  or  Fortran  77. 

Program  POSTCL  calculates  the  return  or  insertion  loss  and  phase 
for  a  given  simple  diameter,  waveguide  size  and  frequency  scan  and 

for  a  given  dielectric  constant  and  conductivity.  Frequency 
dependent  permittivities  can  be  generated  by  other  codes  and  entered 
in  POSTCL  by  a  data  file.  This  program  is  used  to  estimate  the 
experimental  data  from  a  known  sample  (Fig.  4).  For  example,  if  a 

carbon  fiber  sample  of  a  given  diameter  and  known  D.C.  conductivity 

is  to  be  studied,  we  can  calculate  if  there  will  be  enough 

sensitivity  to  detect  it  or  if  it  wi’l  cause  such  a  large  impedance 
mismatch  so  as  to  be  impractical  to  analyze.  One  can  then  pre-select 
samples,  sometimes  bundling  several  lengths  of  material  before  the 
actual  experiments  are  attempted. 

Program  EPOST  takes  return/insertion  loss  and  phase  data  and 
matches  it  by  using  Newton's  method  to  iterate  the  equations  at  each 
freque.rty .  The  program  requires  an  initial  guess  for  the 
permittivity.  One  has  the  option  of  starting  out  the  program  wth  a 
single  guess  for  the  first  frequency  and  allowing  the  prog  >m  to 
precede  to  subsequent  frequencies  starting  from  the  initial  guess  or 
starting  from  the  solution  for  the  previous  frequency,  ^ne  also  has 
the  option  of  specifying  the  initial  guess  at  each  frequency  by  means 
of  ...  input  data  file.  Such  a  file  is  typically  generated  from  an 
analysis  of  a  similar  sample  or  by  extrapolation  from  the  frequency 
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Figure  5.  A  plot  of  the  relative  permittivity  of  water.  The  data 
were  taken  at  X-band  using  the  HP8510  Vector  Network  analyser  and  the 
measurement  was  repeated  three  times  with  the  sample  repositioned 
between  measurements.  The  dots  are  EPOST  inversions  of  the  data  for 
both  Sll  and  S21.  The  crosses  result  from  a  REFPOSTTEK  analysis  on 
the  three  sets  of  data.  The  line  was  generated  by  program  Posttek 
operating  on  all  three  data  sets  and  fitting  the  data  to  second  order 
polynomial  representati ves  for  the  frequency  dependence  of  g'  and  G". 
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dependent  fit  of  the  sample  at  a  different  waveguide  band.  The 
program  generates  plots  (Fig.  5)  and  a  permittivity  file  which  can  be 
used  with  other  programs. 

Program  REFPOSTTEK  performs  a  least  square  fit  to  S-parameter 
data  at  each  frequency  to  find  e.  It  differs  from  EPOST  in  that  a 
least  square  fit  is  used  instead  of  a  Newton  inversion  of  the  data. 
The  least  square  approach  is  less  likely  to  diverge  at  a  bad  initial 
guess  for  e.  Since  both  the  real  and  imaginary  parts  of  e  are  found, 
this  program  requires  at  least  three  data  points  at  each  frequency. 
It  thus  differs  from  EPOST  which  requires  exactly  two  data  points  for 
each  Newton-Rolfson  inversion.  The  two  couples  S-parameters  Sn 
and  S21  result  in  four  data  points  at  each  frequency  and  thus 
suffice  for  the  REFPOSTTEK  analysis.  Often  samples  with  large 
permittivities  and  non-cylndrical  crossections  or  slightly  tapering 
diameters  will  produce  differing  S-parameters  each  time  they  are 
inserted  in  the  waveguide  fixture.  Program  REFPOSTEK  will  “average" 
over  the  scatter  in  measured  data.  REFPOSTEK  is  similar  to  POSTTEK 
(described  below)  in  that  it  performs  a  least  square  of  the  data  but 
it  differs  in  that  the  inversion  is  performed  at  each  frequency.  It 
does  not  require  a  specification  of  a  functional  from  for  e  with 
frequency.  The  input  data  is  taken  from  a  single  waveguide  band 
measurement  series. 

Program  POSTTEK  takes  complex  reflection  and/or  transmission  data 
at  one  or  more  waveguide  bands  and  performs  a  least-squares  best  fit 
of  the  data  to  a  menu  of  selected  complex  functions  of  frequency 
representing  the  frequency  dependence  of  the  permittivity.  The 
program  will  vary  one  or  more  of  the  parameters  appropriate  to  the 
designated  function  and  print  out  the  optimized  values  and  plot  the 
fit  to  the  data  (Fig.  6).  The  optiizaton  of  the  parmeters  is 
performed  by  standard  least  squares  subroutines(S) .  The  menu 
includes  a  variety  of  functions  including  polynomials,  a  conductivity 
function,  a  damped  harmonic  oscillator  function,  a  Lorentzian,  a 
Gaussian,  etc.  The  input  files  are  prepared  from  the  AVNA 
S-parameter  files  or  from  dB  and  phase  files  by  subsidiary  codes 
(NANFIL  and  X YDBPHS2 V ) .  A  dielectric  function  file  may  be  created 
using  the  parameters  and  other  subsidiary  programs  such  as  PERFIL  and 
TRANSXY.  The  permittivity  file  may  serve  for  a  first  guess  for  EPOST 
in  analyzing  new  data  or  for  comparison  to  the  output  of  EPOST. 

Program  FTEMTEK  is  similar  to  POSTTEK  in  that  it  will  fit  a 
multitude  of  data  to  a  specified  function  of  frequency.  This  program 
fits  the  complex  permittivity  calculated  by  EPOST  or  REFPOSTTEK 
rather  than  the  measured  transmission  and  reflection  coefficients.  A 
plot  of  the  fit  generated  by  FTEMTEK  is  shown  in  Fig.  7.  It  is  found 
that  the  POSTTEK  and  FTEMTEK  programs  often  result  in  very  nearly  the 
same  answers  but  the  FTEMTEK  code  requires  much  less  computer  time 
since  it  does  not  have  to  repeatedly  calculate  the  post  impedance 
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Figure  6.  A  plot  generated  by  program  POSTTEK  showing  the  real  and 
imaginary  reflection  (X  and  Ku  bands)  and  the  fit  for  a  conducting 
fiber  sample  (Mi  substituted  phthal ocyani ne  in  Kevlar). 
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equations.  Since  the  input  to  FTEMTEK  is  one  step  removed  from 
actual  measured  data  and  relys  on  the  proper  convergence  of  EPOST  or 
REFPOSTTEK,  POSTTEK  is  the  prefered  code  in  the  case  of  a  difficult 
data  analysis. 

As  mentioned  above,  program  PHASETEK  calculates  a  phase  from  the 
phase  data  acquired  using  the  transmission  bridge.  The  bridge  puts 
out  a  number  of  dB  insertion  loss  files  versus  frequency  for  various 
settings  of  the  calibrated  phase  shifter  in  the  reference  arm  of  the 
bridge  circuit.  (See  Fig.  8a).  PHASETEK  transforms  the  insertion 
loss  data  (dB)  into  scalar  transmission  coefficient  data  and  then 
performs  at  least  squares  fit  of  the  data  to  the  equation  A  + 
B*sin(C)  at  each  frequency  (Fig.  8b).  The  phase  C(f)  are  saved  on  a 
file  and  the  difference  between  the  phase  with  the  sample  present  is 
subtracted  from  the  pnase  without  the  sample  using  program  ASMDXY. 
The  phase  difference  file  is  used  for  the  evaluation  of  the 
dielectric  response. 

Codes  POSTCLTEK,  POSTTEK,  REFPOSTTEK  and  EPOST  rely  on  a 
subroutine  package  POST  which  calculates  the  equivalent  impedances 
and  the  complex  reflection  and  transmission  coefficients.  This 
subroutine  package  is  listed  in  the  Appendix  along  with  the  relevant 
subroutines  and  calling  programs  of  some  of  the  previously  described 
codes. 


EXAMPLES 


We  examine  the  use  of  the  above  programs  in  the  analysis  of 
several  samples.  The  first  sample  examined  is  a  partially 
graphitized  carbon  fiber  designated  DG-129-R.  A  single  fiber  has  a 
diameter  of  the  order  of  one  micrometer  and  produces  a  negligible 

impedance  change  in  the  waveguide.  A  bundle  of  the  order  of  50 
fibers  needed  to  produce  a  reasonable  signal.  We  estimate  the 
sample  diamctui  of  the  fiber  bundle  by  measuring  the  length  and 
weight  of  the  bundle  and  calculating  the  effective  diameter  from  the 
known  density  of  the  material.  The  plots  of  the  return  loss  data 
produced  by  the  AVNA  at  X  and  Ku  bands  are  shown  in  Fig.  9.  The 

periodic  oscillations  are  due  to  the  presence  of  standing  waves  which 
could  not  be  completely  cancelled  or  corrected  for  in  the  error 
correction  software  of  the  model  8409S  AVNA.  The  insertion  loss  and 
the  phase  difference  produced  by  the  transmission  bridge  after  the 
use  of  program  PHASE  and  subtraction  of  the  phase  output  from  the 
empty  holder  are  shown  in  Fig.  10.  The  Ka  band  bridge  data  appear  to 
be  of  higher  quality  tnan  the  AVNA  data  but  are  more  difficult  to 

acquire  than  that  using  the  automated  system.  Standing-wave  damping 
appears  to  be  more  complete. 

In  the  next  figure  (Fig.  11)  we  display  the  output  files  of  EPOST 
for  all  three  waveguide  bands.  The  standing-wave  oscillations  noted 
in  the  AVNA  data  are  transformed  into  oscillations  in  the  apparent 

permittivity  with  respect  to  frequency.  A  clear  trend  is  seen  in  the 
data  consistent  with  an  inverse  frequency  falloff  of  the  imaginary 
part  of  the  dielectric  constant  as  is  expected  from  a  simple 
conductor.  A  smal1  frequency  dependence  is  seen  in  the  real  part  of 
e.  The  solid  lines  in  the  figure  are  the  frequency  dependent 
permittivity  given  by  program  FTEMEK  assuming  a  frequency  independent 
e  coupled  with  a  complex  conductivity:  <j  =  0. 78645-. 00022 j  and  e  = 
2. 886-0. 3242 j.  The  imaginary  part  of  the  conductivity  is  assumed 
responsible  for  the  small  fall-off  in  the  real  part  of  the 
permittivity  with  increasing  frequency,  the  coefficients  calculated 
for  a  similar  functional  fit  to  the  reflection  and  transmission  data 
via  program  POSTTEK  are  very  nearly  the  same  within  a  few  percent. 
Since  the  X  and  Ka  band  data  are  rather  noisy  we  have  also  fit  just 
the  Ka  band  data  by  itself  to  a  frequency  indepedent  complex  e  and  a 
real  conductivity  using  POSTTEK.  This  fit  results  in  the  following 
coefficients:  e  =  13.20  -  187  j  and  a  real  conductivity  a 
=2.560/0hm-cm.  The  conductivity  is  within  12  of  that  inferred  from 
the  larger  set  of  data. 


As  a  second  example  we  study  the  case  of  a  polycrystal ine 
ferroelectric  sample.  A  sheet  of  the  material  was  polished  to  a 
thickness  of  12  mil  and  then  sliced  to  produce  needles  with  a 
rectangular  12  x  12  mil  crosssect'ion.  A  needle  was  selected  and 
measured  at  X  and  Ka  bands  employing  the  Model  3510  AVNA.  Although 
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Figure  9.  A  plot  of  the  return  loss  of  a  conducting  fiber  bundle 
taken  by  the  automated  vector  network  analyser  at  X  and  Ku  bands. 
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Figure  10.  A  plot  of  the  Insertion  loss  and  phase  difference  data 
taken  on  the  transmission  bridge  at  Ka  band  for  conducting  fiber 
bundle  (DG-129-R). 
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Figure  11.  The  EPQST  dielectric  points  along  with  fits  to  the  data 
(lines)  for  fiber  bundle  sample  DG-129-*.  The  data  was  taken  at 
three  different  waveguide  oands;  X,  Ku,  and  Ka.  The  oest  fit  curve 
was  generated  using  parameters  from  FTEM  or  POSTTEK  (nearly  identical 
numbers)  for  all  three  bands. 
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reasonable  results  can  be  obtained  for  sampies  with  small  rectangular 
crossections  where  an  effective  sample  diameter  is  calculated 

assuming  a  circular  crossection  of  equal  area.  No  waveguide 
attenuators  were  used  to  suppress  standing  waves.  The  S  parameter 
data  shows  much  less  scatter  than  encountered  in  the  previous  example 
at  the  same  frequency  bands.  (See  Figure  12).  The  data  was  fit 

using  POSTTEK  and  assuming  a  quadratic  frequency  dependence,  e*  =  A  + 
Bf  +  Cf2)  -  j(0  +  E*  +  Ff2) ,  over  the  frequency  range  studied, 

where  f  is  the  frequency  in  GHz.  A  slightly  different  ferroelectric 
formulation  was  studied  at  X-band  using  the  HP8510  Netcwrk  Analyzer. 
Three  independent  measurements  of  the  X-band  S-parameters  were 
analyzed  using  program  REFPOSTTEK.  The  twelve  input  curves  (Re  and 
Im  Sn  and  Spi )  are  shown  in  Figure  13  along  with  the  inverted 

complex  permittivity. 


CONCLUSIONS 

This  report  describes  the  theory,  experimental  set-ups,  and 
computer  cedes  developed  at  NRL  for  the  use  of  the  dielectric 
post-in-waveguide  technique  for  the  measurement  of  dielectric 
materials.  The  technique  is  especially  useful  in  studying  the 
properties  of  filamentary  high  loss  dielectric  materials.  The  method 
is  mostly  used  for  room  temperature  studies  though  some  work  has  been 
done  warming  the  sample  holder  and  reference  line  to  40C  and  cooling 
to  IOC.  We  expect  to  be  able  to  extend  this  limited  temperature 
range  to  -100  to  +30QC.  Tne  experimental  technique  is  simple;  the 
data  analysis  is  also  strightforeward  but  the  method  requires 
substantial  computer  computations. 
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APPENDIX  I 

Measurement  and  Analysis  Procedure 

In  this  section  we  list  the  steps  involved  in  the  measurment  and 
data  reduction  of  a  sample  by  our  DPW  method. 


SAMPLE  SELECTION  AND  MOUNTING 

A  sample  is  chosen  with  a  large  length  to  diameter  ratio  to 
minimize  depolarization  effects.  Usually  the  sample  is  inserted  in  a 
Debye-Scherer  capillary  tube  with  outside  diameter  small  enough  to 
fit  through  the  smallest  hole  in  the  waveguide  test  fixtures 
(0.5mm).  For  lower  dielectric  constant  material  a  post  is  fashioned 
for  each  waveguide  band  with  diameters  usually  less  than  a  third  of 
the  waveguide  width  and  length  closely  matching  the  waveguide 
height.  The  latter  samples  are  press-fitted  into  position  in  split 
waveguide  sections  instead  of  being  inserted  through  a  hole  in  the 
waveguide.  For  high  conductivity  fibers  or  some  filamentary  linear 
chain  ferroelectric  crystals  a  fine  fiber  filament  is  found  which 
spans  the  waveguide  at  the  lowest  measurement  frequency  and  is 
cemented  across  a  waveguide  cover  flange  with  conducting  paste.  A 
small  grove  is  cut  in  the  flange  o  mount  the  sample  in  this  mannaer 
alligning  it  to  spar-  the  middle  of  the  waveguide.  The  sample  is  cut 
from  the  flange  fixture  with  a  razor  blade  and  mounted  in 
successively  smaller  flange  mounts  for  higher  frequency  measurements. 


AUTOMATED  VECTOR  NETWORK  ANALYSER  MEASUREMENTS 

To  use  the  AVNA  it  must  first  be  calibrated  for  measurements  at  a 
particular  waveguide  band.  The  calibration  procedure  is  described  in 
the  AVNA  manual  and  takes  approximately  10  minutes.  The  calibration 
on  the  HP8409S  AVNA  is  performed  with  5  to  10  dB  isolation  pad 
waveguide  sections  in  place.  No  pads  are  needed  on  the  HP  8510 
AVNA.  The  test  fixture  is  in  place  during  the  transmission  segments 
of  the  calibration  procedure.  We  usually  require  the  return  loss 
from  the  empty  test  fixture  after  calibration  to  be  35  dB  down  or 
more  over  the  frequency  range  before  the  calibration  is  accepted  and 
measurements  begin. 

After  the  calibration,  the  sample  is  placed  in  the  test  fixture 
and  the  S-parameters  are  measured.  The  AVNA  lists  and  plots  the 
S-parameters.  From  the  listing  and  plots  generated  by  the  AVNA  we 
determine  if  the  data  is  adequate  for  further  analysis.  If  not  the 
sample  may  be  repositioned,  a  sample  with  a  larger  smaller 
diameter  may  be  found,  or  other  adjustments  may  be  made  to  lower  the 
standing  wave  interference.  Once  the  data  is  acceptable  it  is  stored 
on  a  floppy  disk  before  other  samples  are  measured.  Later  the 
Hewlett-Packard  floppies  are  transcribed  to  IBM  exchange  format  for 
transport  to  a  VAX  computer  for  analysis. 


TRANSMISSION  BRIDGE  MEASUREMENTS 


The  transmission  bridge  is  assembled  as  shown  in  Fig.  3,  The 
sweeper  is  turned  on  set,  to  sweep  in  free  run  across  the  waveguide 
band  and  allowed  to  stabilize  for  approximately  an  hour  before 
measurements  begir.  Initially  the  reference  arm  is  turned  off  by 
setting  the  attenuator  in  that  arm  to  maximum  and  the  signal  through 
the  sample  arm  is  used  to  normalize  the  scalar  network  analyser.  The 
normalization  is  done  with  an  empty  capillary  tube  in  the  place  of 
the  sample  when  appropriate  and  the  sample  arm  attenuator  is  set  to  7 
dB  to  damp  standing  waves.  The  reference  arm  attenuator  is  then  set 
back  until  the  power  through  the  reference  arm  is  slightly  larger 

than  that  through  the  sample  arm  at  all  frequencies.  The  phase 

shifter  set  to  0  deg;  a  trace  is  acquired  by  the  network  analyser  and 
transfered  to  floppy  disk  on  a  DEC  LSI  computer.  The  phase  shifter 
is  advanced  and  other  traces  are  acquired  until  data  is  taken  for 
phase  settings  from  0  to  340  degrees  in  steps  of  20  degrees. 

After  the  empty  bridge  has  been  characterized  the  sample  is 
inserted  in  the  test  fixture.  The  reference  arm  is  again  turned  off 
and  the  intertion  loss  is  measured  on  the  network  analyzer  and 
transfered  to  disk.  The  reference  arm  is  then  turned  on  and  the 

power  set  so  that  the  reference  level  is  once  again  slightly  larger 
than  the  power  through  the  sample  arm  at  all  frequencies.  The  phase 
shifter  is  set  to  zero  and  traces  are  taken  and  stored  on  disk  at 

successive  phase  shifter  settings  as  done  previously.  Once  the  phase 
data  has  been  acquired,  a  new  sample  is  inserted  and  measured.  Some 
bridge  setups  have  a  lOdB  directional  coupler  in  the  sample  arm  to 
measure  the  scalar  reflection  coefficient  of  the  sample.  To  make 
this  measurement,  a  short  is  inserted  at  the  sample  position  to 
normalize  the  network  anlayzer  and  then  the  short  is  removed  and  the 
sample  is  placed  in  position.  The  scaler  netowork  analyzer  trace  is 
recorded  on  floppy  disc  for  later  analysis. 


COMPUTER  ANALYSIS 

Programs  to  analyze  the  SPW  data  are  run  on  a  VAX  11/780 
computer.  In  the  following  description  of  a  data  reduction  we  assume 

that  a  sample  has  been  measured  using  both  a  Hewlett-Packard 

automated  vector  network  analyser  and  a  transmission  bridge. 

Data  from  the  AVNA  is  brought  to  the  PDP  11/34  on  floppy  disks 
written  in  the  IBM  exchange  fromat.  The  floppy  disk  is  then  carried 
to  the  VAX.  Once  on  the  VAX  the  data  is  placed  into  standard  file 
structures  which  are  callable  by  a  wide  vareity  of  programs.  NANFIL 
takes  the  ASCII  S-parameter  file  and  outputs  dR  versus  frequency  and 

phase  versus  frequency  as  unformatted  binary  "XY"  files.  (The  binary 

files  take  up  less  disk  space.)  Since  the  sample  is  usually 
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positioned  some  distance  behind  the  AVNA  reference  plane  when  the 
measurements  are  performed,  a  program  OFFPHS  is  invoked  to  adjust  the 
phase  files  to  account  for  this  offset.  OB  and  phase  files  are  then 
ready  to  be  run  using  program  EPOST. 

Data  from  the  transmission  bridge  is  written  on  floppy  disks  as 
binary  "XY"  files  in  the  RT -11  format  by  the  LSI -11  computer  at  the 
experimental  station.  The  floppies  are  then  carried  to  the  VAX,  read 
and  reformatted  to  VAX  VMS  files  using  program  RT2VMS.  The  insertion 
loss  data  files  are  directly  usable  by  our  analysis  programs, 
however,  the  phase  data  requires  considerable  manipulation  to  obtain 
the  phase  shift  versus  frequency  required  for  later  analysis.  To 
obtain  the  phase  shift  program  PHASETEK  must  be  invoked  twice;  once 
to  find  the  initial  phase  versus  frequency  for  data  taken  with  an 
empty  sample  fixture,  and  then,  to  obtain  the  pase  versus  frequency 
with  the  sample.  Thus  PHASETEK  is  called  up  and  the  relevant  data 
files  are  read  in  and  an  initial  guess  is  given  for  the  parameters  of 
the  fitting  sine  function.  PHASETEK  will  then  output  phase  vs 
frequency  binary  files.  The  phase  file  for  measurements  with  the 
sample  is  then  subtracted  from  the  file  for  measurements  with  an 
empty  sample  holder  to  obtain  the  relative  phase  shift  versus 
frequency.  The  subtraction  is  done  usinq  program  ASMDXY.  The 
insertion  loss  and  phase  files  are  then  ready  for  EPOST. 

Program  EPOST  calculates  the  permittivity  of  a  sample  from  a  dB 
and  a  phase  value  at  each  frequency.  The  program  internally 
minimizes  the  difference  between  the  measured  real  and  imaginary 
parts  of  the  voltage  reflection  or  transmission  coefficient  and  the 
calculated  value  via  Newton's  method.  OB  and  phase  data  files 
generated  for  a  single  measurement  of  a  sample  in  a  given  waveguide 
are  entered  into  EPOST  along  with  the  sample  diameter  and  waveguide 
size.  The  diameter  is  measured  with  a  micrometer  or  a  calibrated 
traveling  microscope  for  a  single  fiber  sample  and  is  calculated  for 
fiber  bundles  from  a  known  density  and  the  measured  length  and  mass. 
Program  EPOST  requires  a  guess  for  the  permittivity  at  the  first 
frequency  in  order  to  start  off  its  itterative  search  for  the 

permittivity.  If  this  guess  is  not  close  enough  the  interations  may 

diverge.  If  such  divergences  are  a  problem,  the  dB  and  phase  versus 

frequency  may  be  calculated  for  a  number  of  guesses  for  the 

dielectric  constant  and  conductivity  using  program  POSTCL  and 
compared  to  the  measured  data.  Or  it  is  often  more  useful  to  take  a 
subset  of  the  input  data  files  (using  ASSMBL  to  form  the  smaller 
files)  and  run  POSTTEK  with  a  simple  dielectric  function  since  this 
program  is  less  susceptable  to  numerical  divergences.  Once  an 
acceptable  initial  guess  is  found,  EPOST  inverts  the  data  frequency 
by  frequency  and  generates  binary  permittivity  files.  Observation  of 
the  EPOST  plots  will  usually  suggest  trial  function  for  the  frequency 
dependence  of  the  permittivity. 


If  enough  data  is  available  program  REFPOSTTEK  can  be  employed 
instead  of  EPOST.  REFPOSTTEK  performs  a  least  squares  fit  to  the 
S-parameter  data  at  each  frequency  point.  It  requires  more  input 
data  than  fitted  parameters;  thus,  to  calculate  the  real  and 
imaginary  parts  of  the  permittivity  at  least  three  data  curves  are 
required.  They  could  be,  for  instance,  the  real  and  imaginary  parts 
of  S21  along  with  the  magnitude  of  S21 -  Though  more  data  a^e 
prefered  including  several  independent  measurements  of  the  same 
sample  or  several  samples  of  the  same  material  with  different 
diameters.  To  use  REFPOSTTEK  the  dB  and  phase  data  are  rewritten  as 
real  and  imaginary  S-parameters  vs  frequency  via  program  XYDBPHSTOV. 
Amplitude  S-parameter  files  are  created  from  dB  data  by  program 
DB2VXY.  The  least  square  interaction  procedure  is  not  as  sensitive 
to  the  initial  guess  of  the  permittivity  as  is  EPOST.  The  program 
will  accept  ?  permittivity  file  generated  by  another  program  to 
specify  the  initial  values  for  s  at  each  frequency  or  it  will  accept 
an  initial  guess  of  the  real  and  imaginary  parts  of  the  permittivity 
at  the  first  frequency  point  and  enquire  whether  this  guess  is  to  be 
applied  at  all  subsequent  points  or  whether  subsequent  points  are  to 
be  initiated  at  the  previous  soluation. 

Program  POSTTEK  fits  the  measured  complex  reflection  or 
transmission  data  from  one  or  more  measuremnts  of  the  sample  to  a 
menu  selected  function  of  frequency.  In  order  to  invoice  POSTTEK,  the 
dB  and  corrected  phase  files  are  combined  and  written  to  disk  as  real 
and  imaginary  voltage  coefficient  versus  frequency  files  via  program 
XYDBPH32V.  POSTTEK  will  also  accept  scalar  voltage  coefficients  and 
will  simultaneously  analyze  data  taken  at  several  waveguide  bands. 
PHASETEK  optimizes  the  coefficients  of  the  specifiec  functions  and 
can  put  out  files  of  the  calculated  voltage  coefficients  from  the 
optimized  functions.  The  program  runs  rather  slowly  due  to  the 
lengthy  nature  of  the  post  calcuations. 

Often  nearly  identical  functional  dependences  are  obtained  much 
more  rapidly  using  program  FTEMTEK.  To  use  this  program  permittivity 
files  generated  by  EPOST  or  REFPOSTTEK  are  split  into  real  and 
imaginary  parts  using  program  SPLTFL  and  entered  in  FTEMTEK.  FTEMTEK 
then  fits  the  EPOST/ REFPOSTTEK  data  to  the  same  menu  specified 
funtions  of  frequency  and  outputs  files  of  the  calculated  optimized 
fits.  The  outputted  files  can  be  used  by  program  TID  to  plot  the  fit 
and  the  data  in  formats  suitable  for  publication.  The  output  files 
also  can  be  combined  and  used  as  input  files  for  POSTCL  calculations 
or  for  first  guess  permittivity  files  for  further  EPOST  or  REFPOSTTEK 
i tterations. 

The  measurement  and  calculation  procedures  are  thus  seen  to 
require  many  steps,  but  once  the  equipment  and  computer  codes  are  "up 
and  running"  the  analysis  is  essentially  straightforward. 
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APPENDIX  II 
Program  Listings 


INTRODUCTION 

In  this  appendix  are  listed  some  of  the  FORTRAN  programs  used  in 
the  DPW  analysis.  The  first  listing  is  the  subroutine  package 
POST. FOR  which  includes  the  post  impedance  calculations  and  the 
conversion  to  reflection  and  transmission  parameters.  This  is 
followed  by  the  calling  programs  of  POSTCL.FOR  and  EPOST.FOR.  The 
subroutines  under  the  title  FTPOST.FOR  are  presented.  These 
subroutines  are  linked  to  a  general  fitting  program  package 
(FITTEK.FOR)  to  produce  POSTTEK.  Only  a  small  number  of  changes  are 
required  in  FTPOST  to  produce  FTEM.FOR  which  is  linked  to  the  fitting 
routines  to  produce  FTEMTEK.  The  difference  of  files  FTPOST.FOR  and 
FTEM  FOR  is  listed  in  the  last  part  of  this  appendix. 


THE  POST  CALCUATION:  POST. FOR 

The  following  subroutine  package  calculates  the  post  equivalent 
circuit  impendances  and  the  reflection  and  transmission  parameters 
for  a  dielectric  post  centered  in  a  rectangular  waveguide. 
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SUBROUTINE  POST <EZ1. EZ2, A,  D. F,  EFEP) 


CALCULATES  THE  cOUIVALhNT  IMPEDANCES  OF  A  DIELECTRIC  POST 
IN  WAVEGUIDE 


l 

22 


IMPLICIT  COMPLEX  (E) 

EP  »  CONJCIEPEP! 

PIE  *  3.  1415927 

ALPH  *  0. 1047 1 976*F*D 

EBET  *  ALPH*CSOR  T  <  EP ) 

E JOB  -  EB£9J0(EBET) 

EJ1B  *  EBESJ1 (EBET) 

CALL  BESJO( ALPH, JO, AJOA,  V,  I  ERR ) 

CALL  BESJ1 (ALPH, JO, AJ1A, V,  IERR ) 

ALAM  *  30. /F 

ALAMO  *  2*A*ALAM/SQRT( -ALAM*ALAM  +  4*A*A) 

TYPE  *, F, A, D, AJOA, AJ1A, EJOB, EJiB 
A1  *  A/ (2. *ALAMG ) 

A2  *  P IE*D/A 
A2E0  -  A.2*A2 

EQ1  *  Al*( (EJOS/(AJOA#(£BET*AJOA*EJ1B-ALPH*EJOB*AJ1A) ) ) 
1-S0IA, D, ALAM)>ALPH*ALPH*0.  25) 

EQ2  »  4* A 1 *A2SQ/ ( ( ALPH*ALPH*£ J 1 B / ( A J 1 A* ( ALPH*A JOa*£ J 1 B 
1-£BET*EJ0B*AJ1A) ) )-2.  ) 

El  *  CMPLX  <0.  ,  1.  ) 

EZ2  -  -EI*(0.  5*EQ2-»-EGl  ) 

EZ1  *  EI*EQ2 

RETURN 

END 

FUNCTION  SO(A, b, ALAM) 

ABL»2. +A/ALAM 

SO  -  ALOCI 1. 2732395*A/D>-2. 

DO  1  J  *  1, 100 
AJODD  =>  2*J+1 

ACUM  *  < ( 1 .  /SORT <AJODD*AJODD-ABL*ABL> )-l .  /AJODD ) 
CONTINUE 

SO  -  SO  +  2.  *ACUM 

RETURN 

END 


FUNCTION  El (EZ< T,  N) 


USEr  WITH  E3ES j  TO  CALCULATE  COMPLEX  3E3SEL  FC 

IMPLICIT  COMPLEX ;£) 

ST  -  STN(T) 

EST  -  EZ*ST-N*T 
El  =  CMPLX (0.  .  1.  )  . 

EST  ■  EI*£ST 
EE5T  -  CEXP(EST) 

El  -  0.  5*<£EST+1.  /EEST) 

El  -  El/3. 1415P27 

RETURN 

END 

FUNCTION  FKT) 

IMPLICIT  3VTE  td) 

IMPLICIT  COMPLEX  it! 

COMMON  /VARIB/EZ. N. GREAL,  G 
IFC.  NOT.  QREALiCO  TO  1 
FI  »  REAL (El <EZ.  T.  N) ) 

RETURN 

FI  »  A1MAC (El (EZ, 7,  N)  > 

RETURN 

END 


Cj  O  O 


<„ 

c 

c 

c 

c 


FUNCTION  ETPOST i a. D. F. EP) 

CALCULATES  THE  COMPLEX  TRANSMISSION  COEF 
Sl2  FOR  A  POST  IN  A  RECTANGULAR  WAVEGUIDE 

IMPLICIT  COMPLEX  t£> 


CALL  P0ST(EZ1.  EZ2,  A, D, F, EP ) 

D  TYPE  *,  Zl.  Z2\  EZ1 .  EZ2 

C  ETPOST  *  EZ2/(EZ1*EZ1  *  EZ1  +  EZ2  +  2.  *EZI*EZ2> 

ETPOST  *  2.  *EZ2/(EZ1*EZ1  +  2.  * ( EZ1*EZ2+EZ1+EZ2 i +  1 

RETURN 

END 
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FUNCTION  ERPGS7  v  A.  D.  F,  tP  > 


C  CALCULATES  THE  COMPLEX  REFLECTION  COEF 

C 

C  SU  —  FOR  A  DIELECTRIC  POST  IN 

C 

C  RECTANGULAR  WAVEGUIDE 

C 


IMPLICIT  COMPLEX (E ) 

CALL  POST (£21. E22,  A,  D,  F,  EP ) 

El  -  EZ1  ♦  EZ2 
£12  *  E1*E1 
E22  -  EZ2*EZ2 
E3  -  El  +  1. 

E32  -  E3*E3 

ERPOST  *  (E12  -  1.  -  E22 > / (E32  -  E22) 

RETURN 

END 
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j»  a  &>.  ■  j*  >.  r>j 


o  n  c*  n  o  n  r> 


SUBROUTINE  TPGSTD< TD3,  TPH. a, d. f. ep  ) 


CALCULATES  THE  INSERTION  LOSS  AND  PHASE 

PROM  A  DIELECTRIC  POST  IN  MATCHED  RECTANGULAR  UlAUECuIDE 


IMPLICIT  COMPLEX  (£) 

E  *  ETPOSTIA, D, F. EP ) 

D  TYPE  *,  'T'.  E 

TDB  -  -20.  *ALQGl 0 ( CA0S < £ ) ) 

TPH  *  57  2957S*AIMAGiCL0G(E) ) 

RETURN 

END 
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SUBROUTINE  RFQSTDiP.OG.  SPH.  A,  D.  EP  ) 

C 

C  CALCULATES  THE  RETURN  LOSS  AND  PHASE 

C 

C  FROM  A  DIELECTRIC  POST  IN  MATCHED  RECTANGULAR  WAVEGUIDE 

C  -* 


IMPLICIT  COMPLEX  (£) 

E  *  ERPOST (A. D. r , EP ) 

RDB  -  -20.  *ALOG1O(CA0S(E> ) 

RPH  -  57  29579*A IMAG ( CLOG ( £) ) 

RETURN 

END 
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FUNCTION  REPG5T ( £5 ,  A.  D  F,  QT  > 


CALCULATES  THE  REAL  PART  OF  THE 
REFLECTION  OR  TRANSMISSION  COFFFIEIENTS 


IMPLICIT  SVTE  (0) 

IMPLICIT  COMPLEX (E) 

IF ( QT ) CO  TO  1 

REPOST  «  REA<_  ( ERPOST  i  A,  0.  F,  EP  >  ) 
GO  TO  2 

REPOST  *  REAL( ETPOST (A,  0.  F»  EP  >  > 

RETURN 

END 
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J*  "Li*  TLK  /\jN  _M  < V* 


^UNCTION  AIPOST • EF . A, D,  f . QT > 


CALCULATES  THE  IMAGINARY  PART  OF  THE 
REFLECTION  OR  TRANSMISSION  COEFFIEIENTS 


IMPLICIT  BYTE  (Q) 

IMPLICIT  COMPLEX  (E) 

IF  (QT)  00  TO  1 

AIPOST  *  A I MAO (ERPOSTCA.  D.  F.  EP ) ) 
GO  TO  2 

AIPOST  *  AIMAG<ETPOST<A,  0. F.  £P ) ) 

RETURN 

END 


o  o  o  n  r,  r*  o  o  o  r»  rt  o  o  r*  c»  o  «  j  o  o  o  o  ^  o  r,  o  o  r#  o 


FUNCTION  4BPCST  EP , A, D , F, OT ) 


CALCULATES  THE  ABSOLUTE  VALUE  OF  THE 
REFLECTION  AND  TRANSMISSION  COEFFICIENTS 


IMPLICIT  Q  V  f  E ( ' j  > 

IMPLICIT  COMPLEX lE) 

IF  lQT>  CO  TO  1 

ABPQST  =  CABS( ERPOST ' A,  D.  F,  EP ) > 
CO  TO  2 

ABPOST  =  CABSfET?OST(A,  T.  c. EP)  ) 

RETURN 

END 

*-4*4*-*-*  **-44**4*  *4*4--  -4  » 


SUBROUTINE  EBESUO 


COMPLEX  JO  BESSEL  FUNCTION  FROM  A  POLYNOMIAL  APPROXIMATION 

TAKEN  FROM  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS.  . 

US  DEPT  OF  COMMERCE  —  NBS  APPLIED  MATH  SERIEE  55 
EOS  9.  4.  4  AND  9  4  6 


COOD  FOR  RANGE  -j  TO  3  ERROR  LESS  THAN  1.3E-8 
AND  FROM  3  TO  INFINITY  WITH  ERROR  LESS  THAN  4E-9 

IN  PARAENTHESIS  LIKE  PDP  ROUTINE  EXCEPT  JO  AND  E  DO  NOTHING 
SX  IS  COMPLEXES 
SB J  IS  C0MPLEX*8 
IER  IS  ERROR  CONDITION 


subroutine  eebesjocsx,  JO.  SBJ. E.  IER) 

IMPLICIT  DOUBLE  PRECISION  <A,B.T,F.X> 

C  REAL  A< 18) 

C0MPLEX*8  SBJ. SX 

COMPLEX* 16  BJ,  OBJ.  X3B V ,  X,  XBY3,  TH, FO,  THO 
DOUBLE  PRECISION  M13I 

DATA  A/-2.  2499*97D0,  1.  265620SD0.  3163S66D0.  044-47=00 

1, -  0039 a/, 4D0.  OC 021  DO 

2.  00000077DC.  -  005 5274D0. 000095 1 2D0,  00l3723.'D0. 

3-  00072805DO.  00014476D0 

4. -  04166397D0.  -  0000395400.  .  0026357300-  -  0005412500 

5,  -  00029333DG,  00013558D0/ 

X  =  SX 

XBY3  =  X/3. 

IF  (REAL  <  X  > .  LT  -3  ) IER  =2 
IF  ( ACS( X)  CT  3)  CO  TO  1 
CBJ  =DCMPLX(  1.  DO.  0  DO) 

DO  2  J  =  1.  6 
J2  =  2.  *J 

2  BBJ  =  B3J  *  A(  J;*XBy3**J2 

SJ  =  BBJ 
CO  to  3 
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X3SY  *  1  /  XB  v3 

BBJ  =  DCMPLX <0  7O759456D0, 0.  DO) 

TH  =  DCMPlX ( -  7E539916D0. 0.  DO) 

DO  4  J  =  1,6 
UP  =  U  +  6 
JPP  =  J  +  12 

89  J  »  SB  J  ♦  A  (  JP  )*X3BY**U 
4  TH  =  TH  +  A(  JP  P  )  *X3BY**J 

FO  =  BO J 
THO  »  X  +  TH 

QJ  =  F0*CDC0S; THui  'CDSORT(X) 

3  SBJ  *  CMFLX(GJ) 

RETURN 

END 

C  ■4**  *********  "  ^ 44***44««*4*44  4444  »*  -4-J-4-4 

0 

c 

0  SUBROUTINE  BESJ1 

C 

r 

C  COMPLEX  Jl  BESSEL  FUNCTION  FROM  A  POLYNOMIAL  APPROXIMATION 

C 

C  TAKEN  FROM  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS 

C  US  DEPT  OF  COMMERCE  —  NBS  APPLIED  MATH  SERIES  55 

C  EQS  9,  4.  4  AND  9.  4.  6 

C 

C 

C  GOOD  FOR  RANGE  -3  TO  3  ERROR  LESS  THAN  1.3E-3 

C  AND  FROM  3  TO  INFINITY  WITH  ERROR  LESS  THAN  4£-S 

C 

C  IN  PARAENTHESIS  LIKE  PDP  ROUTINE  EXCEPT  Jl  AND  E  DO  NOTHING 

C  SX  IS  C0MPLEX*8 

C  SBJ  IS  CQMPLEX*8 

C  IER  IS  ERROR  CONDITION 

C 


34444414444444444444444444444444»*44444«444444*44«4*4444 I444H 


SUBROUTINE  EEBESJ1 <SX, Jl , SBJ, E, IER ) 

IMPLICIT  DOUBLE  PRECISION  <A,B.T,F,X> 
C0MPLEX*8  SX , S3 J 

COMPLEX* 16  X,  X3Y3,  X3BY, 88 J, BJ, TH, FI, TH1 
REAL  A (20) 

DOUBLE  PRECISION  A (20) 

DATA  A/-.  562499S5D0,  2 1 093573D0, - .  039542S9D0. 


0044331 


1-  00031761D0,  0000 1 109D0. 

2.  00017 1 05D0,  002495 1 1  DO , 
3-2. 356 1 9449D0 , 

4.  124996 1 2D0,  ,  0000565D0, - 
5-  00029 166D0/ 

>:  =  DCMPLX  (  SX  ) 

XBY3  =  X/3. 

IF  !  REAL ( X ) .  LT  -3.  )  IER  =2 
IF  ( CDABS  <  X )  GT.  3'  GO  TO  1 
B8J  =  DCMPLX (0.  DO,  0  DO) 

DO  2  J  =  1,6 
J2  =  2.  *J 

BBJ  =  BBJ  ♦  A( J) »X3V3**J2 
3  J  =  X*  (  .  5  *  3  3  J  ) 

GO  TO  3 

X33Y  =  i.  / XB y3 


79788454>D0  ,  00  0  00  1  54jD0.  01 
001 I3653D0,  -  00020033D0- 


9D0. 

6  59o6  7D0 • 


006378 , 9D0, . 0007434800, 


)7tJ824DG  • 


jv'j  rUVaA  fiLriLW 


-’'LxLiLir.i 


*  '  r-'i 


o  r>  ft  Ci  Ct  o  cj  Cj  o  n  r>  o  o  o  o  o  o  o  o  o  r> 


BBJ  =  DCMPLX ( 0  CO. 0  DO 
TH  =  DCMPLX  (0.  DO.  0.  DO) 

DO  4  J  =  1 , 6 
JP  »  J  *  7 
jpp  >  j  +■  14 

BBJ  =  BBJ  *  A  (  JP  )  •»X3Q  V**J 
4  TH  *  TH  *  A ( JPP ) *X3BY<  *J 

FI  =  A l 7 )  +  OBJ 
TH  1  -  X  +  A  <  1  4  /  *  TH 
BJ  =  FI  -*CDC0S  ( Th 1 ) /CDSGRT ( X ) 

3  SB J  *  BJ 

RETURN 
END 

C  ■*•***  *  -**■*■**-+  *>4444  ■*  *  #  ■#  « IJ-*  *  * 


SUBROUTINE  BESJG 


JO  BESSEL  FUNCTION  FROM  A  POLYNOMIAL  APPROXIMATION 

TAKEN  FROM  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS. 

US  DEPT  OF  COMMERCE  --  NBS  APPLIED  MATH  SERIES  5: 
EOS  9  4  4  AND  9  4.6 


OOOD  FOR  RANGE  -3  TO  3  ERROR  LESS  THAN  1.3E-8 
AND  FROM  3  TO  INFINITY  WITH  ERROR  LESS  THAN  4E-9 

IN  PARAENTHESI5  LIKE  PDP  ROUTINE  EXCEPT  JO  AND  E  DO  NOTHING 
SX  IS  SINGLE  PREC IC I ON  INPUT 
SB J  IS  SINGLE  PRECISION  OUTPUT 
IER  IS  ERROR  CONDITION 


C  8  ■*■*•»■*■**■*  ******  ■*•*■*■***-*■*  ■♦*■*■**•*■*  »■*■*■*•*  ■»■*■*■»*■*-*■»■*■*■*■*-♦  *  ***•»<'* 

SUBROUTINE  BESJO<SX,  JO. SB J. E.  IER  ) 

IMPLICIT  DOUBLE  PRECISION  (A.B.T.F.  X) 

C  real  A( 18) 

DOUBLE  PRECISION  A ( 1 8 ) 

DATA  A/-2.  2499997D0,  1  .  26562U8U0-  -  3163866D0,  044AA7A’D0 
1,  -  0039444D0.  .  00021  DO 

2, 0000O077D0, -  005527400.  -  D0009512D0,  .  00137237P0. 

3-.  00072805D0.  00014476D0 

4,  04 166397D0. -  00003954D 0.  .  OOS62573DO.  -  00054125D0 

5, -  00029333D0.  00013058D0/ 

X  =  SX 

XBY3  *  X/3. 

IF  (  X.  LT.  -3  )  IEP  =2 
IF  ( ABS ( X  ) .  GT  3)  GO  TO  1 
BBJ  =  1. 

DO  2  J  =  1.6 
J2  =  2.  *J 

2  BBJ  =  BBJ  *  A(  Ji*X3Y3-»-»J2 

BJ  =  BBJ 
GO  TO  3 

1  X30Y  =  1  /XBY3 

BB^  =  0  7978845600 
TH  =  -  7853981oDG 
DO  4  j  =  l  ,  6 
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JP  —  J  *  a 
JPP  =  J  -  12 

38 J  =  30 J  ♦  a< JP'*x3cY»>j 
TH  =  TH  -  A(  jPP  ■  *\3BY*-*J 
FO  =  BDJ 
THO  »  X  *■  TH 

8 J  =  FO*DCQS ( THO ) / OSGR  T ( X  ) 

3  SB  J  =  03 

RETURN 
END 

C 

lw 

C  SUBROUTINE  BESJ1 

C 

C 

C  J1  3ESSEL  FUNCTION  FROM  A  POLYNOMIAL  APPROXIMATION 

C 

0  TAKEN  FROM  HANDBOOK  CF  MATHEMATICAL  FUNCTIONS 

C  US  DEPT  OF  COMMERCE  —  NBS  APPLIED  MATH  SERIES  55 

C  EOS  9  4  4  AND  9,  4  & 

C 

c 

C  GOOD  FOR  RANGE  -3  TO  2  ERROR  LESS  THAN  1  3E-8 

C  AND  FROM  3  TO  INFINITY  WITH  ERROP  LESS  THAN  4E-S 

C 

C  IN  PARAENTHES I S  LIKE  POP  ROUTINE  EXCEPT  J1  AND  E  DO  NOTHING 

C  SX  IS  SINGLE  PRECICION  INPUT 

C  SB J  IS  SINGLE  PRECISION  OUTPUT 

C  IER  IS  ERROR  CONDITION 

C 

c 

C3*44*****4*#4*44*4#4«4444**4*#*4*444444 44444444*444 4444444444444444 

SUBROUTINE  BESJt (SX,  Jt, SBJ, E,  IER) 

IMPLICIT  DOUBLE  PRECISION  <A.B,T,F,X) 

C  REAL  A (20) 

DOUBLE  PRECISION  A(20) 

DATA  A/-  56249985DO.  21093572D0,  -  02954289CC.  OC 4A 22  1 9DC , 

1-.  00031761D0.  00001 1 09D0.  7978845oD0,  OOOOOlS&DO.  Olfc59ao~DC 

2.  000 1  7 1 05D0,  -  0024951  IDO,  00U3653D0,  -  0002002300. 

3-2  35a l 9449D0, 

4  12499M2D0.  000056500,  -  00637879D0.  OG074348D0,  00079824DG. 
5-  00029 1 6AD0 / 

X  -  SX 
XBY3  =  X/3. 

IF  <  X.  L7  -3.  )  IE?  =  2 
IF  (ABS< X )  GT  3)  GO  TO  1 
B3J  =  0 
DO  3  U  =  1,  A 
J2  =  2  *J 

2  3BJ  =  3BJ-*  a  i.  j<  J2 

OJ  =  X*<  .  5  +  OBJ) 

GO  T0  3 

1  X30Y  =  l/XBy3 

BBJ  =  0 
TH  =  0 

DO  4  J  =  i , 6 
JP  =  U  »  7 

JPP  =  J  4-  14 

BBJ  =  3BJ  *  x(.P'4x2GV«4j 
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TH  *  TH  +  A(  JPP  >*X3BY**J 
FI  *  A ( 7 )  +  SB  J 
TH  1  *  X  +  A  <  1  4  )  +  TH 
BJ  *  F1*DC0S(TH1)/DSQRT(X> 

3  SB  J  *  BJ 
RETURN 
END 

SUBROUTINE  SIMULA I  I,  VAR,  QK,  GOUT) 
IMPLICIT  BYTE(u) 

DIMENSION  VAR (II) 

DIMENSION  DVaR ( 2 .' ,  EP  ( 2 ) 

GOUT  =  FALSE. 

D  TYPE  *,  'SIMUL ’ 

DO  1  I  *  1 .  1 1 
DVAR ( I )  *  VAR ( I ) *  OOOl 
IF ( DVAR ( I )  EG. 0.  ) DVAR (  I )  =  1.  E-5 
1  CONTINUE 

TYPE  *,  'VAR  *  ' , VAR ,  'DVAR  *',DVAR 
TYPE  *.  'VAR  = ' , VAR,  'DVAR  *'.DVAR 
CALL  CORREC (EP,  1 1,  VAR, DVAR ) 

TYPE  *,  'EP  »•',  EP 
EPSIL  -  0. 

DO  3  I  *  1,  II 
El  «  0. 

IF<EP( I). EQ.  0.  OR  VAR  cl).  EG.  0  >G0  TO  3 
El  *  ABS ( EP ( I ) /VAR (  I  )  ) 

EPSIL  »  EPSIL  +  El 

IF  <DMAX. LT. ABS < VAR < 1>  )  )  GOUT  -  TRUE. 
IF  (EPSIL.LT  OK.  OR.  GOUT)  GO  TO  4 
DO  5  I  *  1,11 
VAR  (  I )  *  VAR  (  I  )  *  EP  (  I  ) 

TYPE  *.  'VAR ( I i  , VAR 
IF  (EPSIL.LT  OH) GO  TO  4 
GO  TO  2 

4  RETURN 
END 


SUBROUTINE  CORRECT?,  II, 
DIMENSION  VAR< II ), DVAfl ( I 
TYPE  *, 'CORREC ■ 

DO  I  1  =  1,2 
EP(  I  >  =  0. 

F<I>  =  FUNK  I,  YAR> 

CALL  FDERW<2,  VAR,  DYAR,  Dl 
DENOM  =  DF(1.  1>+DF(2.  2>- 
IF  (DENOM. EQ.  0)  CO  TO  10 
EP ( 1 )  *  (F (2)*0F ( 1, 2>-F( 
EP(2)  =  (F(  1  >  ■♦DF  (2,  1/-FC 
CONTINUE 

TYPE  *-  -FI.  F2. EPSIL1, EPS 
RETURN 


3UER0UTINE  FutnV!J  ;  .  *,  ^AR.  DVAR,  vr ) 
OIMENSIGN  VAR  (II),  DVAR  i  1 1  > .  DP  (II.  ID 

D  TYPE  *,  'FDERVV' 

DO  2  I  *  l  -  1 1 
DO  1  J  *  1,11 
V  *  VAR(J> 

D  =*  DVAR(J) 

D  TYPE  *.  'FDERVV  V ,  D  =  ' .  V ,  ,  D 

VAR(J)  =  V+0 
F  =  FUN  I  (  I,  VAR  i 

D  TYPE  *-  'F«'.  F 

VAR(J)  -  V  -  D 

_  DF  l  I .  J  >  ■  (F-FUNK  I,  VAR)  >/(2.  *0) 

1  VAR(J)  =  V 

2  CONTINUE 

D  TYPE  *. 'FDERVV  DF ( , )  - • , DF 

RETURN 
END 


y*. 


FUNCTION  FUNK  I.  VAR) 

IMPLICIT  BYTE  (Q) 

IMPLICIT  COMPLEX (E) 

DIMENSION  VAR  ( 2 )  •  "* 

COMMON/VALUE/ VAL< 2), A,  D.  F,  QTRANS.  QTANDR 
TYPE  *, 'FUN I ' 

EP  *  CMPLX ( VAR ( 1 ) . VAR  <  2 ) ) 

IF (QTANDR )CQ  TO  6 
CO  TO  (1,2)1 

FUNI  -  REPOST <£P, A, D. F. QTRANS)  -  VAL < 1  ) 
CO  TO  3 

FUNI  -  AlPOSTfEP. A. D. r.  QTRANS)  -  VAL<2> 
CO  TO  3 
CO  TO  (4,  5)1 

FUNI  *  ABPOSTCEP. A, D.  F,  ,  FALSE.  >-VAL< 1 > 
CO  TO  3 

FUNI  -  AflPOST(EP. A, D. F,  .  TRUE.  >-VAL(2) 

RETURN 

END 


THE  CALCULATED  RESPONSE  FROM  KNOWN  SAMPLE:  POSTCL  FOR 


This  program  listing  calculates  the  reflection  and  transmission 
loss  and  phase  for  a  sample  with  a  known  dielectric  constant  and 
conductivity.  The  permittivity  may  be  specified  over  the  waveguide 
frequency  range  using  a  binary  file  prepaired  by  another  code. 


PROGRAM  POSTCl 
IMPLICIT  BYTE  <G.' 

IMPLICIT  COMPLEX i E ) 

DIMENSION  XC201), Y ( 20 1  ) ,  Z ( 20 1 ) 

CHARACTER*25  QNAM,  ON A ME.  GPLTDS. GNAMD,  GNAMP 


TYPE  11 
PRINT  11 

FORMAT (  '  CALCUALTES  THE  INSERTION  OR  RETURN  LOSS  PHASE  OF 
1'  A  DIELECTRIC  POST  IN  WAVEGUIDE ' ) 


INCLUDE  'PLTDES. FOR ' 


QPLT  *  QUER( 'PLOTS?  ( TEK  4103)', 17) 

IF  (QPLT)CALL  PLOTS 

QTRANS  -  QUER<  'TRANSMISSION  CALCULATIONS  25 > 

PRINT  * 

IF (QTRANS) PR  I  NT  *,  'TRANSMISSION  CALCULATION 
IF  C .  NOT.  QTRANS) PRINT  *,  'REFLECTION  CALCULATION- 
PRINT  * 

TYPE  *.  'WIDTH  OF  GUIDE  (MIL)' 

ACCEPT  *, WIDE 

PRINT  *.  'WIDTH  OF  GUIDE  -'.WIDE. 'MIL¬ 
A-2.  54E-3*WID£ 

TYPE  *  ,  'DIAMETER  OF  POST  (MIL)' 

ACCEPT  *, DIAM 

PRINT  *.  'DIAMETER  OF  POST  -'.DIAM,  'MIL' 

D  =  2.  54E-3*DIAM 

YY  FILE  INS 

r 

QYY  -  QUER('YY  FILE  INPUTS', i5> 

IF ( QYY ) THEN 

TYPE  *.  'ENTER  YY  FILE  NAME?' 

CALL  FY INN ( GNAME, NPT , NPTS,  FL,  FH,  DF.  X,  Y,  Z,  0) 
PRINT  10.QNAME 
GO  TO  13 

END  IF 

FORMAT (  '  INPUT  PERMITTIVITY  FROM  FILE:  ',A23> 


TYPE  #.  'FREQUENCY  RANGE:  FLO,  FH I ,  DILTA  F  t GHZ )  ' 

ACCEPT  *. FL, FH. DF 

PRINT  *,  'FREQUENCIES  FROM ' . FL,  TO'.FH,  'IN  STEPS  OF'.DF,  'GHZ 
NSTEP  -  (FH-FLI/DF  +  : 

IF  (QYY)  GO  TO  12 

TYPE  *,  'DIELECTRIC  CONSTANT  =' 

ACCEPT  *, REP, A IEP 

PRINT  *,  'DIELECTRIC  CONSTANT  -'.PEP.AIEP 
COND  -  0 
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r>  ct  Ci  ^  pj 


GCuND  =  GUER  (  CONDUCT  I VI  TV'-  ;  ; 

IF ( GCQND ) THEN 

TYPE  *,  'CONDUCTIVITY  IN  / QHM-C M  = 

ACCEPT  *, COND 

PRINT  *,  '  COND  ICT I V ITY  =  •  ,  COND.  '/OHM-CM' 

END  IF 

12  PRINT  * 

DO  1  J  *  1, NSTEP 

IF  ( QYY )  THEN 

F  -  X(J) 

EP  *  CMPLX (Y< J) ,  Z<  J)  ) 

ELSE 

F  -  FL  ♦  ( J-l )*DF 

AIMEC  *  AIEP  +  1800.  *CQND/F 

EP  =  CMPLX ( REP . AIMEC) 

X( J)*F 

END  IF 

IF ( QTRANS ) 00  TO  20 

CALL  RPOSTD( Y( J) ,  Z( J> . A, D, F, EP > 

GO  TO  21 

0  CALL  TPOSTD( Y( J) ,  Z( J). A, D. F, EP ) 

1  CONTINUE 

PRINT  *.  J, X( J). Y(J),  2 ( J ) 

PLOT  DB  AND  PHASE 

CALL  PLTSTF ( NSTEP .  X,  Y,  1,  4.  'DI  POST',  '  FREQ.  'IN  LOSS 
C  CALL  PLTSTF(NSTEP,  X,  Z,  1, 0) 

CALL  WNTCPY 

CALL  PLTSTF( NSTEP,  X,  Z,  1,  5,  'DI  POST',  '  FREQ.  ',  '  PHASE 
CALL  WNTCPY 

QSAV  «  GUER (  'SAVE  ON  FILES?'.  14) 

IF (  NOT.  QSAV) GO  TO  2 
PRINT  * 

IF < QTRANS ) TYPE  *,  'INSERTION  LOSS  FILE' 

IF (  NOT. QTRANS) TYPE  *.  'RETURN  LOSS  FILE' 

CALL  X YOUT ( NSTEP ,  QNAMD, A, D, 0.  ,  X, Y, 0) 

PRINT  14,  QNAMD 

14  FORMAT (  '  DB  FILE'.X.A25> 

IF (QTRANS) TYPE  *,  'INSERTION  PHASE  FILE' 

IF (  NOT.  QTRANS) TVPE  *,  'RETURN  PHASE  FILE' 

CALL  XYOUT ( NSTEP ,  QNAMP  ,  A,  D,  0.  ,  X,  Z,  0) 

PRINT  13, QNAMP 

15  FORMAT ('  PHASE  FILE',X,A25) 

2  CONTINUE 

STOP 

END 
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PROGRAM  evaluates  the 

OF  A  POST  INSERTED  IN  THE  CENTER 

OF  HATCHED  RECTANGULAR  WAVEGUIDE 

GIVEN  THE  INSERTION  OR  RETURN  LOSS 

IN  DG  AND  THE  ASSOCIATED  PHASE  SHIFT  IN  DEG 

OR  THE  SCALAR  RETURN  AND  INSERTION  LOSS  IN  DG 


IMPLICIT  BYTE  iG> 

IMPLICIT  COMPLEX  <E ) 

DIMENSION  Ft  501 ) , Y<  301  ).  Z<  301 ) , VAR (2). VARIG'2) 
1, FST (501).  YSTC501 >,  ZST(501 ) 

COMMON  /VALUE/ VAL ( 2 ) . A,  D, FF, GTRANS,  GTANDR 
CHARACTER*25  GDG, GNAME.  GPHASE,  GNAMEI 


CALL  PRDAT 
INPUT  DATA  TYPE 


FORMAT ( '  INSERTION  LOSS  (DG)  XY-FILE  INPUT ’ ) 

FORMAT ('  RETURN  LOSS  ( DB )  XY-FILE  INPUT') 

FORMAT ( '  ASSOCIATED  PHASE  (DEG)  XY-FILE  INPUT') 

FQRMAT( '  THE  FILES  DO  NOT  HAVE  THE  SAME  NUMBER  OF  PTS!  ') 
GTANDR  -  FALSE. 

GTRANS  =»  QUER(  'TRANSMISSION  MEASUREMENTS  ',  25  ) 

IF(.  NOT  GTRANS )  GTANDR  =>  QUER  (  ’SCALAR  T  AND  R  INPUT'-,  21) 


READ  DATA  FILES 


IF (GTANDR) GO  TO  4 
IF (.NOT  GTRANS) GO  TO  7 


TRANSMISSION  DG  AND  PHASE  FILES 


PRINT  ♦,  'PROGRAM  EPOST  — 

1S0LVES  DIEL  CONST  FROM  POST  TRANS  DB  AND  PHASE 
TYPE  11 

CALL  XYIN(NPTS,  GDB,  A,  3,  C, F,  Y, 0) 

TYPE  13 

CALL  XYIN(NPTSS,  GPHASE, A,  3,  C,  F,  Z.  0) 

IF ( NPTS  EG.  NPTSS ) GO  TO  5 
TYPE  14 
GO  TO  8 


REFLECTION  DG  AND  PHASE  FILES 


PRINT*,  'PROGRAM  EPOST  — 

1  SOLVES  DIEL  CONST  FROM  POST  RE^L  DG  AND  PHASE' 
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<j  - 1  n~  tr\  .•  o  o  >.»  «..•  -  oc>  u  u  (_> 


CALL  X  f  IniNPTS,  vDB,  a,  3  C.  F,  Y  O.' 

TYPE  1C 

CALL  XYlNlNPTSS, OPHASE . A, 3, C, F.  Z, G> 

IFfNPTS.  £Q.  NPTSSJGO  TO  5 
TYPE  14 
GO  TO  7 

C  REFLECTION  AND  TRANSMISSION  D3  FILES 

PRINT  *,  PROGRAM  EPOST  — 

1SQLVES  DIEL  CONST  FROM  POST  SCALAR  REFL  AND  TRANS 
TYPE  12 

CALL  XYIN!NPT5,  QDB. A, B, C. F, Y, 0) 

TYPE  1 1 

CALL  X Y I N ( NP TSS, QPHASE, A. B. C, F,  Z,  0) 

IF  (NPTS.  EG.  NPTSSXJO  TO  9 
TYPE  14 
GO  TO  4 
PRINT  > 

PRINT  15, QDB,  QPHASE 

5  FORMAT  ('  iNPUT  FILE  NAMES:  RETURN  LOSS  FILE= ' ,  X . A25,  ' , 

1  INSERTION  LOSS  FI LE- ' ,  X, A25 > 

GO  TO  6 
PRINT  * 

PRINT  10, QDB,  QPHASE 

0  FORMAT ( '  INPUT  FILE  NAMES:  DB  FILE- ' , X, A25, 

1PHASE  FILE- ',  X,  A25) 

PRINT  ♦ 

WAVEGUIDE  AND  SAMPLE  PARAMETERS 

TYPE  *,  WIDTH  OF  WAVEGUIDE' 

ACCEPT  *, WGW 

PRINT  <*.  'WAVEGUIDE  WIDTH-',  WGW,  'MIL' 

TYPE  *,  'DIAMETER  OF  SAMPLE' 

ACCEPT  *, DIAS 

PRINT  *,  'DIAMETER  OF  SAMPLE- ',  D  IAS,  'MIL' 

E  FORMAT!//,  '  STEP  ',  5X ,  'FREQ ',  1  OX,  'FILE  l',7X,  FILE  2'.  8*. 

1  'RE(EP  )  9X,  '  IM(EP  >  ',  //  ) 

CONVERT  TO  CM  FROM  MIL 

D  -  2.  S4E-3*D I AS 
A  -  2.  54E-0*WGW 

INITIAL  GUESS  OF  DIELECTRIC  CONSTANT  AT  LOWEST  FREQ 
QLAST-  FALSE 

QFIL=QUER( 'START  FROM  YY  FILE  INFUTT',25) 

IF(QFIL)THEN 

CALL  FY  INN  <  QNAME I ,  N, NPT,  FI, F2, DF, FST ,  YST ,  ZST, 0) 

PRINT  22, QNAME I 
TYPE  22, QNAME I 
22  •  FORMAT!'  STARTING  VALUES  FROM  -FILE  ',X,A25> 

VAR! 1 ) =YST ( 1 ) 

VAR ( 2 ) =  ZST ( 1 ) 

ELSE 

TYPE  *,  INITIAL  GUESS  OF  DIELECTRIC  CONSTANT  -'T  *  i  1  ,  ' 

ACCEPT  * , VAR 


QLAST  =  QUERl  'START  AT  i_AST  VALVE '  ' ,  20 ) 

IF  t  OLAST ; THEN 

PRINT  *.  'ITERATION  STARTS  FROM  PREVIOUS  FREQ  SOLUTION 
ELSE 

PRINT*,  'ITERATION  STARTS  FROM  INITIAL  GUESS  AT  E.A-Ih  FREQ  ' 

END  IF 
PRINT  IS 

>ARIC(  1  )  »  VAR  i  t  ? 

VARIGT2)  -  VAP  •  > 

ENDIF 

SET  ACCURACY 
OK  =  !  E  -4 

OK  =  (SORT (VAR i 1 i*VAR( 1 > +VAR ( 2 > *VAR ( 2 ) )  >*1.  E- 
El  *  CMPLX  (0  ,  1.  ) 

START  LOOP 

DO  1  J  =  1 >  NPTS 
FF  =«  F  ( J ) 

IF ( OF IL ) THEN 
L=  1 

DDF  -  ABS ( FF-FST ( L ) ) 

IF (DDF. LT  Dr ) THEN 
VAR ( 1 >*YST(L) 

VAR ( 2 ) *ZST ( L ) 

CO  TO  19 
ELSE 

Ic  ( L.  C£.  NPT )  THEN 

TYPE  *,  'NO  FILE  STAKTING  VALUE  FOUND  —  STARTING  FROM  LAST  VALUE 

GO  TO  19 

END  IF 

L»L  +  1 

END  IF 

GO  TO  17 

END  IF 

IF  (QTANDR)CO  TO  20 

EIPH  »  E I *Z ( J ) *1  7453293E-2 

AMAGN  *  10.  *#(-Y( J>/20  > 

EVAL  »  AMAGN#CEXP(£IPH) 

VAL ( 1 )  *  REAL (EVAL ) 

VAL  ( 2  )  *  A IMAG  ( EV'.  .) 

GO  TO  21 

VAL  (  1  )  =  10.  **'.-Y<  J)/20  ) 

VAL ( 2 )  *  10  **<-Z(J)/20  ) 

CALL  SIMU'L(2,  VAR,  OK,  QOUT  ) 

FORMAT ( 15,  5G14  5 ) 

PRINT  40. J.F( J), Y( J),  Z( J), VAR 
TYPE  40.  J, F( J) , Y( J),  Z '  J) . VAR 
Y( J)  =  VAR (  1  ) 

Z(J)  *  VAR (2) 

IF  <  NOT  QLAST  AND  NOT  QF IL) THEN 
VAR ( 1 ) =VAR IG(  1  ) 

VAR ( 2 ) =VAR I G  (  2  > 

END  IF 
CONT INUE 


DF  =  F ' 2 ) -F(  1 

QPL  =  QUER  <  '  PLOTS  7  t  TEk  4105)  ',17) 

IF < QPL ) THEN 

call  plots 

CALL  PLTSTFINPTS. F, Y, l, 1,  'REAL  EP  ' ,  ■  FREQ.  'REAL  EP  '  ) 
CALL  WNTCPY 

CALL  PLTSTFfNPTS, F. Z,  1, 2, '  IM  EP  FREQ  Itt  EP  ') 

CALL  WNTCPY 
END  IF 

TYPE  *.  'DIELECTRIC  CONSTANT  OUTPUT  FILE 
0  TYPE  *,  'REAL  PERMITTIVITY  FILE  OUT ' 

C  CALL  XYOUTfNPTS,  QNAME.  A.  3.  C.  F.  Y,  0) 

C  PRINT  *.  'REAL  PERMITTIVITY  FILE  NAME'. QNAME 

TYPE  *.  'IMAGINARY  PERMITTIVITY  FILE  OUT' 

C  CALL  XYOUTINPTS,  QNAME,  A.  B.  C,  F,  Z,  0) 

C  PRINT  *,  'IMAGINARY  PERMITTIVITY  FILE  NAME '. QNAME 

CALL  FYQUTN( QNAME.  2,  NPTS. F ( l ) ,  F( NPTS ) ,  DF,  F,  Y,  Z) 

TYPE  23, QNAME 
PRINT  23. QNAME 

23  FORMAT ( / ,  '  OUTPUT  FILE  NAME:  ',A25> 

STOP 

END 
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C'JRVF  FITTING  AT  EACH  FREQUENCY:  REFPOSTTEK  FOR 

Program  REFPOSTTEK  Fits,  real,  imaginary  or  magnitude  S-parameter 
post  data  at  each  frequency  in  a  waveguide  band.  The  initial  value 
for  the  fit  is  specified  as  a  complex  constant  employed  at  all 
frequencies  or  employed  at  only  the  initial  frequency  with  subsequent 
starting  values  taken  as  the  solution  at  the  previous  frequency  fit. 
Alternately,  the  initial  guess  can  be  a  permittivity  file  generated 
by  another  program.  Program  REFPOSTTEK  is  less  sensitive  to  be 
initial  guess  than  EPOST  and  will  average  over  multiple  data  sets. 


unnoonn  o  n  o  oooooooooooooooooo 


PROGRAM  REFPostTek 


Least  Square*  fit  of  (Re.  Im,  Abs  vi 1 )  Sll.  S21 
data  on  post  sample  in  waveguide.  Need 
at  laast  three  data  inputs  for  operation. 

For  instance.  th#  input  could  ba 
Ra  and  Ira  Sll  and  S21  data  for  a  post 
sarspla  run  on  tha  HP9310  run  savaral  times 
to  gat  an  avaraga  rasponsa  or  data  from 
savaral  samplas  of  tha  sama  material  at 
different  diameters. 

Input  Files:  A  »  Uavaguida  width 
B  “  Sample  Diameter 
C  »  1  Ra(raf 1  ) 

2  Im<  ref  1  ' 

3  Abs(rafl) 

4  Ra<trans> 

3  Im(trans! 

6  Abs(trans) 

►♦♦♦♦♦♦W************************************* 


BOINQ  IN  DIMENSION.  COMMON.  ETC. 


IMPLICIT  BYTE<G> 

DIMENSION  XINC201 ) >  YIN <201 > . REEP (201 > .  AIMEP (201 > . CHISAV(201 ) 
DIMENSION  CHI <4), XP(2). YP(2> 

DIMENSION  YYDAT (60. 201 > 

COMMON  /FRQ/F(201 ) /HOW/UIDTH( 60 > . DIAM(60 > .  HOU( 60 ) 
1/PAR/PAR/VAR/VAR<2).  DAR<2) 

2/QNAME/QNA<1E(60>.  QNAMFL 
3/LOCAL/YDAT  (60) . YCAL(60> 

CHARACTER+23  QOUT, QPLTDS. QNAM.  QNAME 
DATA  QNAME/ 60+ '  '/ 


TYPE  *.  'PROORAM  RefPostT&k  ' 

PRINT  +.  'PROGRAM  RefPostTek' 
INCLUDE  'PLTDES.  FOR ' 

QPLT*QUER( 'PLOTS?  ( TEK4 103  )  ' ,  1 7 ) 
IF (QPLT)CALL  PLOTS 


gat  input  file  name* 


CALL  NAMES (NFLS.  Al.  A2.  A3) 
QPUNT  =  FALSE. 


c 


Bring  in  files 


DO  1  J  »  l.NFLS 
QNAM-QNAME(J) 

CALL  XYNIN(NPTS, QNAM. UIDTHt J), DIAM( J).  HOW ( J ) ,  XIN.  YIN.  0> 
TYPE  *. J, WIDTH< J).  DIAM( J).  HOW< J) 


IF  <  J.  EQ.  1 )  CO  TO  4 


IF(NPTS.  EQ.  NPT)CO  TO  4 

TYPE  *, 'FILES  HAVE  DIFFERING  NUMBER  OF  POINTS!?' 
CO  TO  3 


NPT  -  NPTS 
DO  3  N-l. NPTS 
IF ( J.  EQ.  1 >F(N>  -  XIN(N) 
DUM  -  YIN(N) 

YYDAT ( J. N)-DUM 


Plot  input  filas 


IF<QPLT) THEN 


IF< J. EQ.  1 ) THEN 

XP(1)  -  X IN( 1 > 

XP (2)  -  X  IN ( NPTS ) 

YP ( 1 )  -  -1 
YP<2)  -  1 

CALL  PLTSTF(2. XP. YP. -1, 3.  'S-Parameter '  Frequency  '.QNAM) 
END  IF 


CALL  PLTSTF(NPTS.  XIN.  YIN.  0.  J.  '  '  '.QNAM) 


END  IF 


CONTINUE 


IF (QPLT)CALL  WNTCPY 


Initial  guess  on  Re  AND  Im  Pts. 


TYPE  *.  'INITIAL  C  JESS  ON  Re(ep)  and  Im(ep)' 

ACCEPT  *. VAR ( I ) ,  V  -  f 2) 

TYPE  *. VAR 

IF  ( VAR  (  l  ).  EQ.  0.  AND.  VAR<2>.  EQ.  0.  )  THEN 

CALL  FYINN < QNAM. N,  NPT, FI, F2,  DF. XIN,  REEP,  AIMEP.  0) 

IF (NPT. NE.  NPTS) THEN 

TYPE  *,  'Fila  numbtr  of  pts  different  from  input  files! 
TYPE  *,  'Try  again ' 

CO  TO  62 

END  IF 

PRINT  61, QNAM 

FORMAT (/.  '  Starting  from  ap  fila;  ' . A25. / ) 


PRINT  *,  'Initial  guess  Re ( ep )  =  ' , VAR ( l ) ,  '  Im ( ep )  =  ' , VAR i 2 ) 
DO  K  >  l, NPTS 


r  VA>\iVa  V>V V.>V£VA--*>\  s  ’ 


END  IF 


END  DO 


KEEP (K ) *VAR  < 1) 

A I MEP ( K ) »VAR  ( 2  ) 


TYPE  *.  'TYPE  incr«m*nts' 

ACCEPT  *. DAP  !  1  > , DAR ( 2 ) 

TYPE  #.  DAR 
NTERMS  -  2 

QLAST  -  QUER < 'Start  from  last  valu*?',22> 
IF!QLAST)PRINT  *.  'Starting  from  previous  » 
PRINT  * 

C 

QPLP-.  FALSE. 

IF!QPLT)QPLP-QU£R<  'PLOT  IN  LOOP?  ' ,  13 ) 

C 

C  START  FREQUENCY  /  FIT  LOOP 

C 

NFREE  -  NFLS  -  NTERMS 
C 

C  SET  UP  FREQUENCY  LOOP 
C 

NPT  -  0 
JFS  -  1 
NFSTEP  -  1 
JFFF  -  NPTS 


C  FREQ  OFFSET  AND  STEP 

C  IF  < I SENSE < & > )  18.  18.  19 

19  TYPE  *.  'STARTING  FREQUENCY  *« ' 

ACCEPT  22.  JFS 

TYPE  *.  '  STARTING  AT  FREQUENCY ',  F( JFS) 
22  FORMAT! IS) 

TYPE  *.  '  FREQUENCY  STEP  INCREMENT  =» ' 
ACCEPT  27. NFSTEP 
27  FORMAT! 18) 

TYPE  *, 'INCREMENT  »', NFSTEP 

C  CLEAR  CHI  SQUARE  ACCUMULATOR 

18  DO  14  J  *  1.4 

16  CHIIJ)  *  0. 

C  BEGIN  FREQUENCY  LOOP 

DO  40  JF  *  JFS, NPTS, NFSTEP 

NPT  -  NPT  1 

IF!. NOT.  QLAST) THEN 

VAR !  1  )  *  REEPUF) 

VAR ! 2 )  -  AIMEP(JF) 

END  IF 

C  DOWN  LOAD  DATA  AT  FRQ  JF 


1  1 
C 


DO  11  K  «  1, NFLS 
DUM  ■  YYDAT ! K, JF ) 
YDAT(K )*DUM 


mm 


w . 
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V 


l'jtiion.  ' 


C  SET  UP  FIT  LOOP 

C 

ALAMDA  -  00 1 

MODE  -  0. 

FRQ  -  F(JF> 

C 

C  PLOT  DATA  IN  LOOP 

C 

IF(. NOT.  OPLP >00  TO  20 
IF ( ( (NPT -1 ) /3) *3+1 . EQ.  NP T ) THEN 
NPC  -  7 

CALL  PLTSTFINFLS, HOW, YDAT,  1. -NPC>  '  '  How  eodt  '  U  ') 


NPC  -  NPT  -  (NPT/S)*3  +  6 

CALL  PLTSTF (NFLS.  HOW.  YDAT,  0,. -NPC.  '  '  ') 

END  IF 
C 

C  TOP  OF  FIT  LOOP 

C 

20  CONTINUE 

CALL  RPOSTCUR (FRQ< ALAMDA,  CHI SQR.  NFLS,  NTERMS,  0) 

C  SEE  IF  CHI  SQUARE  IS  CHANOINO 

CHH  -  0. 

DO  31  K  -  2. 4 

CHH  -  CHH  ♦  CHI(K) 

KKK  -  K  -  1 

31  CHI (KKK)  -  CHI(K) 

CHI (4)  *  CHI SQR 
CHH  -  CHH  ♦  CHI SQR 

IF  (CHISQR.  LT.  1.  E-30)  GO  TO  30 
CHH  -  A8S(CHH/(4. *CHISQR>-1.  ) 

30  CONTINUE 

32  IF  ( I SENSE (2)  EQ. 0)  GO  TO  33 

TYPE  *,  'CHISQR CHISQR,  'LAMBDA', ALAMDA,  'R«(*p)  and  Im(#p)'.VAR 
IF  ( ISENSE ( 4 ) .  EQ.  0 )  CO  TO  33 

TYPE  34,  (YDAT (K) , YCAL (K) >  K  «  l.NFLS) 

34  FORMAT  ( 4(2013.  4, 4X ) ) 

33  CONTINUE 

C  TEST  CHI  SQUARE 

IF  (  (CHH.  LE.  1.  E-3).  OR.  (CHI  (4).  LE.  1.  E-6).  OR.  (ALAMDA.  EQ.  1.  2E6)  ) 

1  GO  TO  13 
GO  TO  20 

C  END  FIT  LOOP 

13  CONTINUE 

C  PLOT  AFTER  FIT 

IF  < .  NOT. QPLP ) GO  TO  35 
NPC  »  NPC  -  6 

CALL  PLTSTF (NFLS, HOW, YCAL.  0, NPC,  '  '  '  > 


*7  ,y[ 


7  „  x 

«■  »  *  J  * 


T  >0  K.”  -  ",  V  *  "  O  4  * *  *  -  "  -  *  *  -  - , 

'Vu’VV'C.  ■ •  ■»'</* <*»  .\V/“V 


35  TYPE  #, JF, CHISQR,  F<JF), VAR 

C  SAVE  Re  and  Im  tp  and  CHISQUARE 

REEP(NPT)  -  VAR ( 1 > 

AIMEP(NPT)  -  VAR ( 2 ) 

CHISAV(NPT)  -  CHISQR 

C  END  FREQUENCY  LOOP 

40  CONTINUE 


C  OUTPUT 


50  CONTINUE 
CALL  TYDAT 
CALL  PR DAT 

PRINT  *,  'PROGRAM  REFPostTek  OUTPUT i 
PRINT  5 1 i QNAMFL 

51  FORMAT (  '  INPUT  NAME  FILE  -'.X,A25> 

PRINT  * 

PRINT  52,  <QNAM£<K>, K»l,  NFLS) 

52  FORMAT ( 5 ( X, A25 ) ) 

PRINT  53 

53  FORMAT (5X,  '#',  10X,  'CHI-SQUARE',  4X,  'Freq.  ',  10X,  'Re  Cap  )  9X,  'Im(ep  )  '  ) 

DO  34  J  *  l.NPTS 

54  PRINT  *, J, CHISAV ( J ) ,  XIN( J )  ■  REEP ( J) ,  AIMEP ( J ) 

IF (QPLT) THEN 

CALL  PLTSTF<NPT,  XIN. REEP,  1,3,  'Re <ep>  ','  FREQ.  ','Re<ep>  ') 
CALL  WNTCPY 

CALL  PLTSTFJNPT.  XIN,  AIMEP,  1,  1.  'Im<ep  FREQ.  Im(*p)  ' 
CALL  WNTCPY 

END  IF 

c 

PRINT  * 

PRINT  *,  'Permittivity  output  file' 

FI  -  F(JFS) 

DF  -  F ( JF5+NFSTEP > -F 1 

CALL  FYOUTN(QOUT,  2,  NPT,  FI,  FRQ,  DF,  F, REEP, AIMEP) 

TYPE  53, QOUT 
PRINT  55- QOUT 

53  FORMAT (  '  OUTPUT  FILE  NAME  -'.2X.A25. X) 

CALL  CODE( 2 ) 

STOP 

END 


FUNCTION  RPOSTFCT< XI, I ) 

COMPLEX  EP 
LOGICAL  Q 
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COMMON  /VAR/VAR (2), DAS (2) 
1/HQW/WIDTH<  60 ) , DI AM< 60 ) . HOW ( 60 > 

EP  *  CMPLX ( VAR  < 1 ) . VAR ( 2 ) ) 

W  -  WIDTH ( I )*2.  34E-3 
D  -  DIAM(I )*2.  34E-3 
Q  -  FALSE. 

IH  -  HOW (  I  ) 

00  TO  (1.  2,  3.  4,  3.  6)  IH 

4  0  ■  .  TRUE. 

1  RPQSTFCT  -  REPOST < EP , U,  D. XI.  Q) 

00  TO  10 

3  0  -  .  TRUE. 

2  RP QSTFCT  -  AIPOSTIEP, W, D,  XI, Q> 

00  TO  10 

6  Q  -  .  TRUE. 

3  RPOSTFCT  -  ABPOST(EP, W,  D, XI,  Q> 

10  RETURN 

END 


SUBROUTINE  RPOSTCUR<FRQ, ALAMDA,  CHirQR.  NFLS.  NTERMS. MODE) 

REAL* 16  ARRAY, WEIGHT, DET. ASQ. ALPHA.  BETA 
COMMON  /VAR /VAR (2) .  DAR  <  2 ) 

2/LQQP/ARRAY(2«  2) 

3/DER/DER IV(2  J 
3/LOCAL/YDAT  <  60 ) . YCAL (60 ) 

DIMENSION  BETA<2>.  ALPHA <2. 2) . YI (60) 

DIMENSION  A(2). B(2),  C(2) 


C  INCLUDE  'EMUL.  FOR' 


C  IF  < I SENSE ( 13 ) .  EQ.  1 )  TYPE  140 

C  140  FORMAT (  '  C  '  ) 

CALL  ERRSNS 

D  TYPE  *,  'CUR  IN', ALAMDA, CHISQR,  NFLS,  NTERMS,  MODE 


l: 


11  NFREE  -  NFLS  -  NTERMS 


61  DO  62  I  -  l.NFLS 
YCAL1I)  -  RPOSTFCT !FRQ,  I  ) 

62  YI(I)  -  YOAT(I)  -  YCAL(X) 

63  CHISQ1  -  FCHISR1NFLS. NFREE,  MODE) 

31  DO  34  J-l.NTERMS 
BETA1J)  -  0. 

DO  34  K-l.J 
34  ALPHA! J,  K)  -  0. 


ii.% 


41  DO  SO  I-  l.NFLS 

CALL  FDER(FRQ, I. NTERMS) 

DO  46  J- 1 .  NTERMS 

BETA!  J )  -  BETA!  J)  +  YI  !  I  >*DERIV!  J  ) 
D  TYPE  *.  J.  BETA! J).  YI! I),  DERIV(J) 


DO  46  K-l. J 

ALPHA!  J.  K)  -  ALPHA!  J.  K)  DER I V ( J )  *DER  I V !  K  ) 

46  CONTINUE 

50  CONTINUE 

51  DO  53  J-l.NTERMS 
DO  S3  K-l,  J 


S3 


ALPHA !K. J)  -  ALPHA! J. K) 


71  DO  74  J-l, NTERMS 

DO  ^3  K-l,  J 

ASQ  "  SORT ! ALPHA ! J,  J ) *ALPHA ! K,  K ) ) 

ARRAY !  J,  K )  -  ALPHA! J, K> /ASQ 

73  ARRAY ! K, J >  -  ARRAY! J. K) 

74  ARRAY !  J.  J )  - 1 .  + ALAMDA 
80  CALL  MAT IRR(NFLS. NTERMS.  DET) 

81  DO  86  J-l.NTERMS 

A! J)  -  VAR ! J ) 

B  !  J )  -  A !  J ) 

DO  86  K- 1 , NTERMS 

96  B!J)  -  B!J)+BETA(K)*ARRAY!J,  K ) /SORT ' ALPHA  <  J,  J ) * ALPHA  v  K .  K )  ) 


47 

D 


DO  12C  J  =  1,  NTERMS 
C(J>  =  A ! J ) 
TYPE  *,  B  <  J) 


61 


VAR  <  J  >  -  A  (  J  ) 

120  CONTINUE 


91  DO  92  I-l.NFLS 

YCAL<I)  -  RPOSTFCTCFRQ, I > 

92  CONTINUE 

93  CH ISQR  -  FCHISR<NFLSi NFREE.  MODE) 

CALL  ERRTST < 10<  ITST) 

IF  (CHISOl  -  CHISQR)  93,  101.  101 

93  ALAMDA  -  10,  *ALAMDA 

DO  130  J  -  1.  NTERMS 

A< J)  -  C(J) 

VAR  <  J )  -  A  <  J ) 

130  CONTINUE 

ICURCY  -  ICURCY  *  1 
IF  <  ICURCY.  OE.  80  >  00  TO  102 
IF  < ITST. EQ.  1)  00  TO  102 
IF  (ALAMDA.  OE.  1E19J00  TO  102 
GO  TO  71 

102  CHISQR  -  CHISQ1 

ALAMDA  *  1 .  2E7 

101  CONTINUE 

ALAMDA  »  ALAMDA/ 10. 

110  RETURN 
END 


SUBROUTINE  FDER(FRQ.  I.  NTERMS) 


C 


C  SUBROUTINE  FDERIV  (NON  ANALYTICAL) 

C  PURPOSE 

C  EVALUATE  DERIVATIVES  OF  FUNCTION  FOR  LEAST-SQUARES  SEARCH 

C  FOR  ARBITRARY  FUNCITGN  OIVEN  BY  FUNCTN 

C  USAGE 

C  CALL  FDERIV  ( X,  I , A, DAR , NTERMS,  DER I V > 

C  DESCRIPTION  OF  PARAMETERS 

C  X  -  ARRAY  OF  DATA  POINTS  FOR  INDEPENDENT  VARIABLE 


oooot^oo  on  oo*~  o  no  nan  noooo 


I  -  INDEX  OF  DATA  POINTS 

A  -  ARRAY  OF  PARAMETERS 

OAA  -  ARRAY  CF  PARAMETER  INCREMENTS 
NTERMS  -  NUMBER  OF  PARAMETERS 
DERIV  -  DERIVATIVES  OF  FUNCTION 

SUBROUTINES  AND  FUNCTION  SUBPROBRAMS  REQUIRED 
FUNCTN  (X,  I,  A) 

EVALUATES  THE  FITTING  FUNCTION  FOR  THE  I  TH  TERM 


MODIFIED  FOR  NET  ANALYSIS  --  SEE  DERIVE  &  YFTE'S 


u-**u-*u-*  •*■**» 


COMMON  /VAR/VAR<2) ,  DAR (2) 
1/DER/DERIV<2) 


IF( I SENSEI  13).  EG.  1 )  TYPE  1 
FORMAT ( 'O' ) 


11  DO  18  J  -  1,  NTERMS 
AJ  -  VAR(J) 

DELTA  -  DAR(U) 

VAR(J)  -  AU  +  DELTA 
YFT  -  RPOSTFCT (FRO.  I ) 


VAR(J)  -  AU  -  DELTA 

DERIV(J)  -  (YFT  -  RPOSTFCT  (FRQ,  I ',  >  /  (2.  *DELTA ) 


IF  (DERIV(U) .  EQ-  0.  )  DERIV(J)  *  . OOOCOl 
VAR  ( J  >  »  AU 


22  CONTINUE 

RETURN 
END 

C  FUNCTION  FCHISR 
C 

C  PURPOSE 

C  EVALUATE  REDUCED  CHI  SQUARE  FOR  FIT  TO  DATA 

C  FCHISR  =»  SUM< (Y-YCAL)**2)/SIGMA#*2'/NFREE 

C 

C  USAGE 

C  RESULT  *  FCHISR  < Y.  S IGMAY, NFLS, NFREE, MODE,  YCAL > 
C 

FUNCTION  FCHISR (NFLS, NFREE,  MODE) 

DOUBLE  PRECISION  CHISQ 
COMMON  /LOCAL /YD AT ( 60 ) , YCAL ( 60 ) 

IF  ( I SENSE! 13)  EQ.  I ) TYPE  1 


1  FO,  MAT  (  '  X  '  ' 

11  CHISQ  -  0 

12  IF  (NFREE)  13,  13,  20 

13  FCHISR  -  0. 

CO  TO  40 

C 

C  ACCUMULATE  CHI  SQUARE 
C 

20  CONTINUE 

DO  30  I  -  i.NFLS 
C  21  IF  (MODE)  22,27,22 

C  22  IF  <  A8S ( YDAT ( I ) ) .  LE.  1 . E-6 >  00  TO  27 
C  IF  (YDAT ( I  )  )  23,27,23 

C  23  WEIGHT  -  1.  /YDAT (I) 

C  CO  TO  30 

C  23  WEIGHT  -  1.  / ( -YDAT ( I ) ) 

C  CO  TO  30 

C  27  WEIGHT  -  1. 

C  CO  TO  30 

C  29  WEIGHT  -  1,  /SIGMAYDAT( I ) **2 

D  TYPE  *.  'CHI  ',  CHISQ, WEIGHT ( I > , YDAT ( I ) , YCAL( I ) 

-30  CHISQ  -  CHISQ  +  ( YDAT( I ) -YCAL ( I > >**2 

C 

C  DIVIDE  BY  NUMBER  OF  DEGREES  OF  FREEDOM 
C 

31  FREE  -  NFREE 

32  FCHISR  -  CHISQ/ (FREE) 

40  CONTINUE 

D  TYPE2, NFREE,  FCHISR 

D2  FORMAT ( 18, Cl3.  3) 

RETURN 

END 


SUBROUTINE  MAT IRR ( NPTS,  NTERMS, DET! 
CQMMON/LQQP/ARRAY ( 2,  2) 

REAL*1&  ARRAY,  AMAX, SA'  .  DET 
DIMENSION  IK(2 ) , JK  ( 2  < 

NORDER  -  NTERMS 
IF  ( I SENSE ( 1 3 ) .  EQ.  1 )  TYPE  1 
1  FORMAT ( *  M' ) 

10  DET  -  1. 

11  DO  100  K-  1, NORDER 
AMAX  -  0. 

21  DO  30  I  -  K,  NORDER 
DO  30  J-K,  NORDER 

23  IF  ( ABS( AMAX > -ABS( ARRAY ( I,  J)  )  )  24,24,30 

24  AMAX  -  ARRAY ( I,  J) 

IK (K ) ■  I 
JK(K)-J 

30  CONTINUE 

31  IF  (AMAX  )  41.32,41 

32  DET  -  0. 

GO  TO  140 

41  I  -  IK (K ) 

IF  (  I-t'. )  21,  51,  43 

43  DO  30  J-l, NORDER 
SAV  *  ARRAY  <K.  J) 

ARRAY  (K, J )  -  ARRAY (  I ,  J  > 


64 


50  ARRAY ( I, J)  -  -SAV 

51  J  -  JK (K ) 

IF  (J-K)  21,61.53 
53  DO  60  t“l. NORDER 
SAV  -  ARRAY! I,  K) 

ARRAY!  I,  K)  -  ARRAY ! I .  J' 

60  ARRAY! I, J)  -  -SAV 

61  DO  70  1=1, NORDER 
IF  (I-K)  63,70,63 

63  ARRAY! I, K)  -  -ARRAY! I , K ) /AMAX 

70  CONTINUE 

71  DO  SO  t-1,  NORDER 
DO  90  J-  1, NORDER 
IF  (I-K)  74,90,74 

74  IF  (J-K)  75,90,75 

75  ARRAY! I, J)  -  ARRAY ( I • J >  +  ARRAY ( I , K> *ARRAY (K, J ) 
30  CONTINUE 

81  DO  90  J»l. NORDER 
IF  (J-K)  93,90,83 
83  ARRAY(K, J)  -  ARRAY (K, J ) /AMAX 
90  CONTINUE 

ARRAY ( K, K )  *  l./AHAX 
D  T'  PE  #,  '  DET'.CET,  '  AMAX  ' ,  AMAX 

100  DET  *  DET*ANAX 

101  DO  130  L-l. NORDER 
K  *  NORDER  -L  +  1 
J  -  IK (K ) 

IF  (J-K)  111,  111,  105 

ios  no  no  :«i.  norder 

?..<V  -  ARRAY!  I,  K) 

ARRAY (  I ,  K )  'i  -ARRAY!  I,  J) 

110  ARRAY! I, 6)  -  SAV 

111  I  «  JK !K ) 

IF  ! I-K)  130, 130, 113 
113  ’JO  120  J  -  1.  NORDER 
•JAV  -  ARRAY (K,  J) 

ARRAY (K,  J )  -  -ARRAY! I,  J) 

120  ARRAY! I. J)  -  SAV 
130  CONYTNUE 
140  RETUHN 
END 


FITTING  MULTIPLE  DATA  SETS  TO  A  FUNCTION  Cf  FREQUENCY:  FTPOST.FOR 

The  following  subroutines  are  linked  to  a  general  fitting  program 
to  produce  POCTTEK. 


bb 


*  “  r  '  Tk  w  »  fwVrf rjC  'TmJf'U' fk,  WV  * *'  -  *f%*  <-  ^  V*.  - 


r  V"  S 

*  ~  m  *  •  *  »  “  %  r  m  14  * 


SUBROUTINE  FTOPT 


r'TDDOP  FOR 


IMPLICIT  QVTE(<?) 

COMMON/OPTION/ I AFCT< 5),  IBFCT<5>, NRT. NCV<5> .  GIND,  OOPT 

1/OPT2/QEPC,  QQQ 

1/NUM/NPTS.  NFLSi  NTERMSi  NPAR 

3/DBV/QDBTOV 

DIMENSION  IFC T ( 25 ) 

CHARACTER*50  FCTLAB 

DATA  IFCT/3.  3,  0,  3,  0,  3.  0,  3.  3,  5,  0,  5.  0,  I,  3.  0,  4,  3,  3,  7,  0,  5,  0 


QOPT  -  .  TRUE. 

QEPC  -  . FALSE. 

GD'JTOV  -  .  FALSE. 

QDBTOV , ■  QUER< 'CONVERT  DB  DATA  TO  V?',27> 


CALL  FTANG(NFLS,  NRT) 


DO  5  J  »  1,23 

CALL  FTFLAB ( FCTLAB.  J> 

TYPE  300, FCTLAB 


NTOT  -  0 
IBFCT(l)  -  1 

DO  I  J  >  1,2 
IF  (J.  EG.  DTHEN 

TYPE  *,  'REAL  PERMITTIVITY' 

PRINT  'REAL  PERMITTIVITY' 

ELSE 

TYPE  *,  'IMAGINARY  PERMITTIVITY' 
PRINT  *,  'IMAGINARY  PFRMI TTIV I TY ' 

END  IF 

TYPE  *.  'FUNCTION  NUMBER' 

ACCEPT  2,  IAFCT(J) 

IA  *  IAFCT(J) 

FORMAT ( IS) 

TYPE  *■  'FUNCTION  #',  IAFCT(U)  - 
CALL  FTFLAB(FCTLAB,  IAFCT(J)) 

TYPE  300, FCTLAB 

PRINT  *,  'FUNCTION  #'.IAFCT<u> 

PRINT  300, FCTLAB 
FORMAT (  X, A60) 


,  0,  0/ 


IF  (  IA.  NE  6.  AND.  IA.  NE  12.  AND.  IA  NE.  15  AND.  Ia  NE  22)  GO  TO  2 
TYPE  *. 'NUMBER  OF  CURVES' 

PRINT  *.  -NUMBER  OF  CURVES' 

ACCEPT  2,  NCV(U) 

TYPE  *,NCV(J> 

PRINT  *, NCV ( J ) 

QEPC  -  QUERC'ADD  A  COMPLEX  CONST  EP^',23) 

C  QIND(IA)  -  QUER (  ’  INDEPENDENT  AMPLITUDES? ',  23 > 

3  I~( IA.  EQ.  10. OR.  IA.  EQ.  12) TYPE*.  'AMPL. 4PIM. X.  LAM.  GAM- 

IF ( IA.  EG.  22)TYPE  *.  'HEIGHT.  ROLL  OFF  FRO.  RELAX  FRO,  ALPHA,  BETA' 
IF <  IA.  EQ.  2.  OR.  IA,  EQ.  4.  OR.  IA.  EQ.  6.  OR.  IA,  EQ.  15>TYPE  *, 

1 'HEIGHT,  CTR  FRQ,  WIDTH' 

NVARA  *  NCV< J) *IFCT( IA) 

IF  (QEPC.  AND.  (  IA.  EQ  6  OR.  IA.  EQ.  12.  OR.  IA.  EQ.  15.  OR.  IA.  EQ.  22)  ) 

1 NVARA  -  NVARA  +  2 
TYPE  *.  NVARA,  'P  ARAME'i  ERS  ' 

NTOT  -  NTOT  +  NVARA 
IF  (J. EQ. 4)  QO  TO  4 
J1  ■  J  +  1 

IBFCT(Jl)  *  IBFCT < J )  *  NVARA 

4  IF ( NVARA. NE.  0)  00  TO  1 

JM1  *  J  -  1 

IBFCT ( J )  *  IBFCT(JMl) 

NCV(J)  »  NCV ( JM1 ) 

C  Q IND ( J )  -  QINDI JM1 ) 

1  CONTINUE 

TYPE  *,  'TOTAL  PARAMETERS  «'.NYOT 
NPAR  *  NTOT 

D  TYPE  *,  '  IAFCT,  IBFCT.  NRT,  NCV.  QIND  '.  IAFCT,  IBFCT,  NRT ,  NCV,  QIND 

RETURN 
END 

SUBROUTINE  FTFLAB ( QFCT,  J ) 

CHARACTER* 30  QFCTN 
CHARACTER* ( * )  QFCT 
QFCTN* ' 

QFCT-QFCTN 

GO  TO  (101,  102.  103.  104,  105,  106,  107,  103,  109,  110,  111,  112,  113,  114,  U5 
1,  1.16.  117,  113.  119,  120,  121.  122.  123)J 
101-  QFCT  n  '  FUNCTN  1:  A  ♦  B*FRQ  ♦  C*FRQ**2' 

GQ  TO  1000 

102  QFCT  *  '  FUNCTN  2:  REAL  LORENTZIAN  ♦  1' 

GO  TO  1000 

103  QFCT  *  '  FUNCTN  3:  I  MAG  LORENTZIAN' 

GO  TO  1000 

104  QFCT  -  '  FUNCTN  4  REAL  GAUSSIAN' 

GO  TO  1000 

103  QFCT  »  '  FUNCTN  3:  IMAGINARY  GAUSSIAN' 

GO  TO  1000 

106  QFCT  *  '  FUNCTN  6:  SEVERAL  REAL  GAUSSIANS' 

GO  TO  1000 

107—  QFCT  -  '  FUNCTN  7:  SEVERAL  IMAGINARY  GAUSSIANS' 

GQ  TO  1000 

103  QFCT  -  '  FUNCTN  8:  CONDUCTIVE  LOSS' 

GQ  TO  1000 

109  QFCT  *  '  FUNCTN  9:  SPOT  DIELECTRIC' 

GO  TO  1000 


nnn  n  n  non 


•.  :0 
1 1  1 
:  12 

1 13 

114 

115 

116 
1  17 
1  IS 

119 

120 

121 

122 

123 

1000 


OFCT  =  FUNCTN  10:  SEAL  KITTEL  PERMEABILITY 
GO  TO  1000 

OFCT  =  '  FUNCTN  11  IMAGINARY  KITTEL  PERMEAO  I L  I  T  v 

GO  TO  1000 

OFCT  -  '  FUNCTN  12:  SEVERAL  REAL  KITTEL' 

GO  TO  1000 

OFCT  -  '  FUNCTN  13:  SEVERAL  IMAGINARY  KITTEL' 

GO  TO  1000 

OFCT  -  '  FUNCTN  14:  CONSTANT' 

GO  TO  1000 

OFCT  -  '  FUNCTN  15:  SEVERAL  REAL  LQRENTZ IANS ' 

GO  TO  1000 

OFCT  -  '  FUNCTN  16:  SEVERAL  IMAGINARY  LQRENTZ IANS ' 

GO  TO  1000 

OFCT  -  '  FUNCTN  17  A  +  3* < FRQ-D)  +  C* ( FRQ-D ) **2 ' 

GO  TO  1000 

OFCT  -  '  FUNCTN  IS.  A  ♦  3*FRQ#*C ' 

GO  TO  1000 

OFCT  -  '  FUNCTN  19:  A  ♦  3*EXP (FRG*C ) ' 

GO  TO  1000 

OFCT  -  '  FUNCTN  20:  PAR ( 1 )  >  RE  DAMPED  HAR  OSC ( 3. 4,  5, 6. 7 >  ' 
GO  TO  1000 

OFCT  -  '  FUNCTN  21.  PAR (2)  IMAG  DAM  HAR  OSC' 

GO  TO  1000 

OFCT  -  '  FUNCTN  22:  SEV  RE  DAMP  HAR  OSC ( H.  f o,  f r , a ,  b  )  ' 

GO  TO  1000 

OFCT  «  '  FUNCTN  23:  SEV  IM  DAMP  HAR  OSC' 

RETURN 

END 


SUBROUTINE  FTQPFN ( I . N, VALU > 

IMPLICIT  9YTE  <Q> 

C  IMPLICIT  DOUBLE  PRECISION  (U) 

IMPLICIT  COMPLEX  <£> 

COMPLEX  KITMU,  DAMHAR 

COMMON/ I NDEP/XI 1001 >, Y< 1001 > ,  YFIT( 1001 ) 

1 /PAR/PAR (51)  t 

2/OPTION/  IAFCT  ( 5  > ,  IBFCT(5>,  NRT  ,  NCV<5>,  QIND,  ODUM 
3/0PT2/QEPC,  QQQ 
DIMENSION  CON( 12 ) 


NTM  -  ISFCT <  N) 

I GO  -  IAFCT (N) 

NCRV  -  NCV(N) 

QINDD  -  QIND' ICO ) 

TYPE  *.  ICO 
NTM2  »  NTM  +  1 

NTM3  -  NTM  +  2 

NTM4  »  NTM  +  3 

NTM5  -  NTM  +  4 

NTM6  -  NTM  +  5 

NTM7  -  NTM  +6 

XI  -  X  <  I  ) 

GO  TO  (  1. 2.  3,  4.  5,  6.  6.  7,  3.  9,  10.  6.  6,  50.  6,  6.  1 1.  12-  13.  14,  15.  6,  4> )  I  GO 
CONST 


o  o  n  *-  r>  r>  n  r>  o  o  ~  n  o  o  -o  n  ci  o  nuinnr.  ^000  n  u  o  o  o  fi  o  (0000  o  *-  «->  o  n  n  r,  •» 


Q  VALU  =  P  AR  1  NTil ) 

IF  (QCYC )  NTM  =  NTM  *  1 
GO  TO  100 

POLYNOMIAL 

1  XI  *  XI  -  PAR ( NTM4 ) 

VALU  *  PAR(NTM)  +  PAR(NTM2)*XI  +  PAR (NTM3) *X  I*X  I 
Ic  (QCYC)  NTM  »  NTM  +  3 
GO  TO  100 

REAL  LORENTZIAN 

VALU  -  1.  +  ALOR  <  X 1/  PAR (NTM) .  PAR  <NTM2) .  PAR  ( NTM3) ,  TRUE.) 

GO  TO  100 

NO  ADVANCE  NTM 

IMAGINARY  LORENTZIAN 

VALU  -  ALOR  <  X I <  PAR (NTM ) . PAR ( NTM2 ) . PAR ( NTM3 ) , .  FALSE.  ) 

IF  (QCYC)  NTM  *  NTM  +  3 
CO  TO  100 

REAL  GAUSSIAN 

VALU  *  1.  *  COMCAU (XI. P AR  <  NTM ) . PAR ( NTM2) . PAR ( NTM3 )  . .  TRUE.  ) 

GO  TO  100 

IMAGINARY  GAUSSIAN 

VALU  -  COMCAU <  XI*  PAR  < NTM)  1  PAR ( NTM2 ) .  PAR ( NTM3 ) ,  .  FALSE.  > 

IF  (QCYC)  NTM  -  NTM  +  3 
GO  TO  100 

MTTEL  PERMEABILITY 

VALU  -  1  +  PAR (NTM)* 

1REAL (KI TMU (XI.  PAR ( NTM2 ) > P AR ( NTM3 > . PAR (NTM4 ) , PAR ( NTM5 )  ) ) 

GO  TO  IOO 

IMAGINARY  KITTEL  PERMEABILITY 
0  VALU  *  PAR (NTM ) * 

1  1A  IMAG  ( KITMU  (XI.  PAR  (NTM2  ) .  PAR  (  NTM3  >  >  PAR(NTM4),  PAR  (  NTM5  )  )  ) 

IF  (QCYC)  NTM  -  NTM  +  4 
CO  TO  100 

DAMPED  HARMONIC  OSCILLATOR 
REAL 

4  VALU  »  PAR ( NTM )  +  REAL ( DAMHAR ( X  I ,  PAR ( NTM3 ) ,  P AR ( NTM4 ) , P AR ( NTM5 ) 

iP AR ( NTM6 ) .  PAR ( NTM7 ) ) ) 

GO  TO  100 

IMAGINARY 

15  VALU  -  P  AR  (NTM2) +AIMAG  (DAMHAR  (  XI.  PAR  (NTM3).  PAR  (NTi-14),  PAR  (NTM5). 

1PAR (NTM6 ) .  PAR ( NTM7 ) )  ) 


u  n  n  n  wooo  a^oo^ooo  ^oon  cdooo  n  m  o  o  o  O'O 


NTM  a 


C  IF  ;GCtL)  NTn  = 

SO  TO  100 

MULTIPLE  GAUSS  IANS 

CAL,.  FITMUL< XI,  VALU,  IGU.  NTM,  NCRV  > 

GO  TO  100 

CONDUCTIVE  LOSS 

VALU  »  PAR ( NTM )  +  1 800.  *PAR < NTM2 > / X I  +■  PAR ( NTM3 )  /  (  X  I *SQRT (XI)) 
IF  (QCVC )  NTM  -  NTM  +■  3 
GO  TO  100 

SPOT  DIELECTRIC  -  RUFUS  WRIGHT 

VALU  »  4.  *PAR  <  NTM ) / <  3.  14 1 5927*PAR < NTM2 > *P AR < NTM3 > > 

GO  TO  100 

POUIER  LAW 

VALU  »  PAR ( NTM )  +  PAR  <  NTM2 )  *  (  X I  **P AR  (' NTM3 )  ) 

GO  TO  100 

EXPONENTIAL 

VALU  ■  PAR  <NTM)  >  P AR < NTM2 > *EXP ( XI #PAR < NTM3 > > 


CONTINUE 

TYPE  *,  XI,  VALU 

RETURN 

END 


FUNCTION  ALOR ( XI , A, B<  C, GRE IM> 
IMPLICIT  BYTE  <Q> 

IF  (. NOT  QREIM)  GO  TO  3 


REAL  LORENTZ IAN 

XX  -  B  -  XI 
D  =*  1.  /C 

ALOR  =  A*XX*D/( 1+XX*XX*D*D> 
GO  TO  100 

IMAGINARY  LORENTZ IAN 


XX  »  B  -  XI 
D  -  1/C 

ALOR  »  A/( 1+XX*XX*D*D) 
100'  RETURN 
END 
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FUNCTION  CQMGau l  X  J ,  A.  3,  C#  ORE ) 

IMPLICIT  3YTE  >  0 ) 

IMPLICIT  DOUBLE  PRECISION  (U) 

CD  TYPE  *,  XI,  A,  B,  C.  QRE 

IF  ( ORE )  GO  TO  4 
GO  TO  5 
C 

C  REAL  GAUSSIAN 
C 

4  XX  *  XI  -  B 
SIGRT2  -  0.  600S#C 
XXBYS  -  XX/SICRT2 

IF < ABS < XXBYS ) . CT  3)  GO  TO  402 
XXBYS2  *  XXBYS#XXBYS 
UXBYS2  -  XX8YS2 
SA  -  -1.  12838*A»XX3YS 
UDSB  -  DEXP < -UXBYS2 ) 

ITM  -  XXBYS2  +  4 
UDSC  -  1. 

DUN  -  1. 

ACLUXX  -  1. 

DO  401  U*l.  ITM 
DUN  -  DUN  #  U 
ACLUXX  -  ACLUXX*XXBYS2 
SD  *  2*U  +  1 

401  UDSC  -  UDSC  +  ACLUXX/ <DUN*SD) 

COMGAU  -  UDSB*UDSC 

COMCAU  -  S.4*C0MGAU 

COMGAU  -  COMCAU  ♦  0.  7979*A*S IGRT2/ ( X I  ♦  B) 
GO  TO  100 

402  COMCAU  -  -1.  12838*A*B/( (B  +  XI)*XXBYS) 

GO  TO  100 

C 

C  IMAGINARY  GAUSSIAN 
C 

5  XX  -  XI  -  B 
SIGRT2  -  0.  6*C 
XXBYS  -  XX/SICRT2 
XXBYS2  -  XXBYS*XX8YS 
COMCAU  -  A*EXP  < -XXBYS2 ) 

IF  (QCYC-)  NTM  =  NTM  3 
CO  TO  100 

C 

C  * 

C 

100  CONTINUE 

CD  TYPE  *. COMGAU 

RETURN 
END 


FUNCTION  DAMHAR (F,  H,  FO,  FR.  DIS,  BETA) 

COMPLEX  DAMHAR, El 
El  =  CMPLX  ( 0  ,  1.  ) 

A  -  H*FR*FO 

DAMHAR  «  A/  ( FO*FO  —  F*F  -  ,EI*FR*F)**(  1.  -DIS)  )-»-*5ETA 
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C*  ■  a  ""a  C  a  _  *  -  **.  a  -  I 

L*  *a>.  j*  ■  "  >> s. 


,V 


»  " 


* » <  s-  r 


' J  <J  <J  O  uug 


RETURN 

END 


FUNCTION  KITMU(F,  Al,  K,  LAM,  GAM) 

REAL  K.  LAM 

COMPLEX  KITMU,  FACl, FAC 2 
A5-Al/<4.  #3.  1415927) 

T 1*2.  #K/A5 

T2-- F*LAM/(A5*GAM#*2) 

FAC  1*CMPLX  ( ( T 1  i-A  1  ) ,  T2) 

FAC2-CMPLX(T1, T2> 

K I TMU-A 1 *FAC 1 /  <  FAc2*FAC 1 - ( F/GAM  >  **2 ) 

RETURN 

END 


FUNCTION  KITM(F,  Al,  AK,  ALAM,  GAM, GR I  > 

REAL  AND  IMAGINARY  KITTEL  FCTS 
OIL  KITMU 

IMPLICIT  8YTE( Q) 

REAL  KITM 
COMPLEX  KITMU 

IF  (QRI)  KITM  *  REAL(KITMU(F,  Al, AK, ALAM, GAM) ) 

IF  (  NOT.  QRI)  KITM  *  AIMAG (KITMU ( F, Al , AK. ALAM,  GAM )  ) 

TYPE  *,  'K',  KITM,  QRI 

RETURN 

END 


SUBROUTINE  FITMUL ( X I ,  VALU,  I GO,  NTM, NCRV) 

IMPLICIT  COMPLEX  (E) 

IMPLICIT  BYTE  (Q) 

COMMON  /PAR/PAR ( 51 ) 

1/QPT2/QEPC, QQQ 
COMPLEX  DAMHAR 

MULTIPLE  GAUSSIANS,  LORENTZ IANS,  KITTEL  FCTS,  DHO  S 

QLOC  -  .  FALjE 
QK IT  -  FALSE. 

QCAU  *  FALSE. 

QDHO  -  .  FALSE. 

IF  ( ICO.  EQ.  6.  OR.  ICO.  EQ.  7)  GO  TO  b 
IF  (IG0.EQ.10.0R.IG0.EQ.il)  GO  TO  40 
IF  ( ICO.  EQ.  15.  OR  ICO.  EQ.  16)  GO  TO  11 
IF  (IGO.  EQ.  22.  OR.  IGO.  EQ.  23)  GO  TO  41 
GO  TO  100 


4i  3DHG  =  TSUE 

KPAR  »  5 
GO  TO  401 

40  GK  IT  -  TRUE. 

KPAR  ■  5 
GO  TO  401 

11  QLOR  -  .  TRUE. 

GO  TO  402 

6  GGAU  *  .  TRUE. 

402  KPAR  -  3 

401  ASM  -  0. 

SSM  -  0. 

OTP  -  .  TRUE. 

IF  ( I  GO.  EQ.  7.  OR.  IGQ.  EG.  11.  OR.  IGO.  EQ.  16.  QR.  IGO.  EG.  23)  QIP  =*  FALSE 


36 


403 


31 


DO  31  J  *  1. NCRU 
JJ  *  J  -  1 

NT  EM  -  NTM  KPAR*JJ 
NTEM2  *  NTEM  +  1 
NTEM3  *  NTEM2  1 
NTEM4  *  NTEM3  +  1 
NTEM3  »  NTEM 4  +  1 
NTEM6  -  NTEM 5  +  1 
AHT  -  PAR (NTEM) 

IF  (QGAU)BSM  -  BSM  ♦  COMGAU( X  I , AHT. PAR (NTEM2 )  ,  PAR ( NTEM3 >  ,  GTP ) 
IF  (QLOR)  BSM  -  BSM  *  ALOR  (  X I .  AHT.  PAR  <  NTEM2)  ,  PAR  ( NTEM3 ) .  GTP  t 
IF  (QKIT) 

1BSM-BSM+AHT*KITM(XI,  PAR (NTEM2 ) . PAR (NTEM3) .  PAR (NTEM4 ) 

2, PAR(NTEM5), GTP) 

IF  (QDHQ ) THEN 

EDHO  *  DAMHAR (XI.  AHT, PAR (NTEM2) ,  PAR ( NTEM3 ) .  PAR (NTEM4 ) 
1,  PAR (NTEM3) ) 

IF (GTP  >THEN 

DHO  -  REAL (EDHO) 


ELSE 

DHO  *  A I  MAG  (EDHO) 

END  IF 

BSM  *  BSM  ♦  DHO 

END  IF 
CONTINUE 


32  IF  (GTP)  BSM  -  BSM  +  1 . 

4 

IF  (QEPC)THEN 

IF (GTP ) THEN 

BSM  »  BSM  -  l.  +  PAR(NTM+KPAR*NCRV) 

ELSE 

.  BSM  =  BSM  PAR<NTM+KPAR*NCRV+1 ) 

END  IF 

END  IF 


VAI..U  <*  BSM 

100  RETURN 

END 

SUBROUTINE  FTANG ( NFLS.  NRT ) 
1  RETURN 

END 
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/>  V  >  w  r, 


SUBROUTINE  FTLAB  < GP ) 

BYTE  GP 

IF ( QP ) PR INT  *,  POST  IN  WAVEGUIDE  ANALYSIS' 
TYPE  *,  'POST  IN  WAVEGUIDE  ANALYSIS' 

RETURN 

END 

FUNCTION  FUN< I ) 


IMPLICIT  COMPLEX  <E) 

IMPLICIT  BYTE  <Q) 

COMMON  / INDEP / X ( 1001 ),  Y< 1001 ).  YFITI 1001 ) 

1 /CONTRL/ AA ( 5 l ) i BB < 51 > . CC < 51 > 

2/NUM/NPTS.  NFLS. NTERMS, NPARV 
3/F^PT/NPTC51  >,  NPTTSS<  1001) 

4/PAR/PAR ( 5 1 ) 

5/OPTION/ IAFCT(5),  IBFCT(5> , NRT, NCV(5l ), QIND(5l ),  GCPT 
6/OFFSET/DD (51) 

XI  *  X ( I  ) 

WIDTH  -  AA ( NPTTSS ( I ) ) 

DIAM  -  BB(NPTTSS( I )  ) 

I  TYPE  =*  CC  (NPTTSS<  I )  )  ♦  .01 


CALL  FTOPFN< I,  1,  A) 
CALL  FTOPFN< I, 2, fl) 
EP  -  CMPLX (A. B) 
CALL  FTQPFN( I. 3, A) 
CALL  FTOPFN< I, 4. B) 
EMU  *  CMPLX (A, B) 


A  »  WIDTH*2.  54E-3 
D  -  DIAM*2.  34E-3 


Q  *  . FALSE. 

GO  TO  (11,  12,  13,  14,  15,  16 ) I  TYPE 
Q  *  .  TRUE. 

FUN  *  REPOSTfEP, A, D, XI, Q) 

CO  TO  20 
Q  -  .  TRUE. 

FUN  -  AIPOSTfEP. A,  D.  XI,  0> 

GO  TO  20 
Q  -  TRUE. 

FUN  *  ABPOST(EP, A,  D,  XI, Q) 

CONTINUE 

RETURN 
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Fi la  SYSPUSERDISKB: CRACHF0RD3FTEM.  FOR; 47 

21  OOBTOV  -  OUER<  ‘CONVERT  OB  DATA  TO  VOt  TACE?  * .  27 > 

22  C 

****** 

Fi  ia  SYSPUSERD I SKB : CRACHFQRD3FTP0ST.  FOR;  24 

21  0D8T0V  -  '1U£R<  ‘CONVERT  OB  DATA  TO  V?'.27> 

22  C 

********** ** 


File  SYSPU3ERDISXB:  CRACHFQRDIFTEM.  FOR;  47 

408  300  FORMAT (  '  FUNCTIONAL  FIT  TO  EP  OR  MU  DATA') 

489  IF(OP JPRINT  300 

490  TYPE  300 

491  RETURN 


MM«« 

Fi la  SYSPUSERDISKB: CRACHFORD 3FTP0ST.  FOR;  24 

488  IF ( OP ) PRINT  a.  'POST  IN  UAVEOUIDE  ANALYSIS' 

489  TYPE  ♦.  'POST  IN  WAVEGUIDE  ANALYSIS' 

490  RETURN 


Fi la  SYSPUSERD I SKB : CRACHF0RD3FTEM.  FOR;  47 


314 

IF  ( I TYPE.  EQ.  1  >  THEN 

313 

CALl.  FTOPFNl  I,  l.  A) 

316 

FUN  -  A 

317 

ELSE 

318 

CALL  FTCPFNd.  2.  B> 

319 

FUN  •  B 

320 

END  IF 

321 

RETURN 

322 

END 

323 

SUBROUTINE  FTSPOUT 

****** 

Fi la  SYSPUSERDISKB: CRACHFORD 3  TPOST.  FOR;  24 


313 

1 

CAIL  FTOPFNd.  1.  A) 

314 

CALL  FT0PFN(I,2.  B> 

313 

EP  -  CMPLX  <A.  B ) 

316 

C 

CALL  FTOPFNU.  3.  A) 

317 

C 

CALL  FT0PFN<I,4.B> 

318 

c 

EMU  -  CMPLX (A.  B )  ^ 

319 

c 

320 

c 

321 

c 

322 

c 

323 

c 

324 

A  -  WIDTH* 2.  34E-3 

323 

0  -  DIAM*2.  34E-3 

326 

c 

327 

c 

328 

c 

329 

0  -  .  FALSE. 

330 

OO  TO  <11.  12,  13,  14,  13,  16 ) ITYPE 

331 

14 

a  -  .  TRUE. 

332 

11 

FUN  -  REPOSTtEP.  A,  0.  XI. 0) 

333 

90  TO  20 

334 

13 

a  -  .  TRUE. 

333 

12 

FUN  -  AIPOSTCEP. A,  D. XI,  0> 

336 

00  TO  20 

337 

16 

a  -  .  TRUE. 
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338 

13 

FUN  -  ABPOSTtEP. A, D. XI. Q) 

339 

20 

CONTINUE 

340 

RETURN 

341 

END 

342 

C 

SUBROUTINE  FT SPOUT 

343 

C 

RETURN 

344 

c 

END 

343 

SUBROUTINE  FTSPOUT 

Nuribar  of  difforonco  faction*  found:  3 
Nuab or  of  difforonco  rocord*  found:  3d> 

DIFFERENCES  / IGNORE- <) /MERCED- 1 /OUTPUT H_P:  CRACHF0RD3FTER  PJFi- 
SYSfUSERDISKB:  CRACHFORDJFTEN.  FOR;  47- 
SYS4USER DISKS:  CRACHFQRDIFTPOST.  FOR:  24 
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